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J.  N.  Reddy 

Department  of  Mechanical  Engineering 
Texas  A&M  University 
College  Station,  Texas  77843-3123,  USA 


Preface 

This  final  technical  report  contains  a  summary  of  research  carried  out  on  Grant  DAAH04- 
96-1-0080  from  the  Army  Research  Office  (ARO)  to  Texas  A&M  University,  with  Dr.  J.  N. 
Reddy  as  the  principal  investigator.  Following  Grantee  data  on  the  project  on  pages  ii- 
V,  the  remainder  of  the  report  is  devoted  to  the  technical  discussion.  To  limit  the  size 
of  the  report  to  a  reasonable  number  of  pages,  certain  topics  are  not  covered  in  detail. 
Additional  information  on  constitutive  modeling  of  shape  memory  alloys,  electrostrictive  and 
magnetostrictive  materials  and  associated  finite  element  formulations  can  be  found  in  the 
dissertation  of  Dr.  Govind  Rengarajan  {On  the  Inelastic  Behavior  of  Crystalline  Solids^ 
Ph.D.  Thesis,  by  Govind  Rengarajan,  Department  of  Mechanical  Engineering,  Texas  A&M 
University,  September  1998). 
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Technical  Discussion 


1.  Background 

In  the  last  two  decades,  the  subject  area  of  smart/intelligent  materials  and  structures 
has  experienced  tremendous  growth  in  terms  of  research  and  development.  Numerous 
conferences,  workshops,  and  journals  dedicated  to  smart  materials  and  structures  stand 
testimony  to  this  growth.  The  vastness  of  the  literature  is  obvious  considering  the 
interdisciplinary  nature  of  the  subject.  Physicists,  mathematicians,  and  engineers  from 
aerospace,  chemical,  civil,  electrical,  materials,  and  mechanical  engineering  fields  are  all 
involved  in  some  part  of  the  development  of  smart  materials  and  structural  systems.  One 
reason  for  this  activity  is  that  it  may  be  possible  to  create  certain  types  of  structures 
and  systems  capable  of  adapting  to  or  correcting  for  changing  operating  conditions.  The 
advantage  of  incorporating  these  special  types  of  materials  into  the  structure  is  that  the 
sensing  and  actuating  mechanism  becomes  part  of  the  structme  by  sensing  and  actuating 
strains  directly.  These  types  of  mechanisms  are  referred  to  as  strain  sensing  and  actuating 
(SSA). 

The  technological  implications  of  this  class  of  materials  are  immense:  structures  that 
monitor  their  own  health,  particularly  useful  in  remote  operations,  process  monitoring, 
vibration  isolation  and  control,  and  medical  applications,  to  name  only  a  few.  On  the 
threshold  of  21st  century,  we  are  entering  into  research  and  development  of  the  next  generation 
of  smart  materials  and  structural  systems.  The  next  generation  of  smart  material  systems  will 
feature  thermo-electro-mechanical  coupling,  functionality,  intelligence,  and  miniaturization 
(down  to  nano  length  scales).  With  the  advent  of  these  new  generation  of  materials,  reliability 
and  integrity  of  these  systems  become  central  issues.  These  systems  operate  under  varying 
conditions  -  they  span  the  whole  spectrum  of  magneto-electro-thermomechanical  conditions. 
These  conditions  could  vary  from  low  to  high  temperatures,  low  to  high  pressures,  low  to 
high  load  levels,  low  to  high  strain  levels,  and  low  to  high  electric  and  magnetic  fields. 
Such  operating  environments  pose  serious  problems  to  the  design  and  maintenance  of  the 
smart  structural  systems.  Experimental  investigations  of  both  the  smart  material  and  the 
structural  system,  though  possible,  are  prohibitively  expensive,  and  therefore  they  should 
be  complemented  with  theoretical  analyses.  The  present  study  is  concerned  with  one  such 
analysis.  The  following  review  of  literature  provides  a  background  for  the  study. 

2.  Review  of  literature 

2.1  Smart  materials  and  structures 

The  phrase  smart  structural  system  refers  to  a  wide  variety  of  active  material  and  passive 
structural  systems.  For  instance,  a  sufficiently  general  system  is  a  composite  (beam,  plate, 
shell,  or  any  other  fundamental  form)  with  embedded  or  surface  mounted  piezoelectric  or 
electrostrictive  patches,  or  even  layers  of  active  materials  in  a  laminated  system  (see  Figure 
2.1).  The  literature  contains  many  definitions  of  smart  and  structures.  We  prefer  to 

follow  Newnham’s  definitions^.  The  structures  with  surface  mounted  or  embedded  sensors  and 
actuators,  that  have  the  capability  to  sense  and  take  corrective  action  are  referred  to  as  smart 
structures.  The  feedback  circuitry  linking  sensing  and  actuating  is  external  to  the  sensor  and 
actuator  components.  In  fact,  this  precisely  distinguishes  a  smart  structural  system  from 
an  intelligent  structural  system.  Intelligent  structural  systems  involve  smart  components  in 
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which  the  functions  of  sensing,  feedback  control,  and  actuating  are  all  integrated.  This  type 
of  system  finds  applications  in  aircraft  wings,  helicopter  rotors,  weapon  systems,  automobiles, 
and  so  on. 


Figure  2.1:  A  schematic  of  a  laminated 


plate  with  surface  mounted  piezoelectric  patches. 


Essentially,  a  smart  structural  system  is  a  multi-functional  unit  (see  Figure  2.2).  It  has  a 
load  (electrical,  thermal,  magnetic,  or  mechanical)  bearing  part  which  is  usually  passive,  and 
an  active  material  part  that  performs  the  operations  of  sensing  and  actuating.  For  example, 
vibration  amplitudes  in  a  flexible  plate  structure  may  be  suppressed  using  the  sensing  and 
actuation  capabilities  of  piezoelectric  or  piezoceramic  films  by  bonding  them  to  the  surfaces 
of  the  plate.  As  the  plate  deforms  due  to  external  applied  loads,  the  bonded  piezoelectric 
film  (sensor)  also  deforms,  and  due  to  its  constitutive  behavior,  it  develops  a  surface  charge 
proportional  to  the  applied  force.  The  charge  may  be  processed  by  a  control  system,  which 
supplies  an  appropriate  voltage  to  the  piezolectric  film  (actuator)  that  induces  a  counteractive 
deformation  to  the  plate  structure  and  the  amplitudes  of  vibrations  are  suppressed. 


Figure  2.2:  Components  of  a  smart  structural  system. 
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Piezoelectric  materials,  electrostrictive  materials,  magnetostrictive  materials,  shape 
memory  alloys,  and  electro-rheological  fluids  are  some  of  the  smart  materials  available  today. 
Among  the  currently  available  sensors  and  actuators,  the  smallest  ones  are  of  the  order 
of  few  millimeters.  However,  progress  towards  intelligent  structures  requires  us  to  develop 
smart  material  systems  that  are  of  the  order  of  a  few  microns.  The  reduction  in  size  has 
tremendous  technological  benefits;  however,  clear  understanding  of  reliability  and  system 
integrity  are  vital  to  the  efficient  and  optimum  use  of  these  material  systems.  As  dimensions 
get  smaller,  induced  electro-thermo-mechanical  fields  get  larger.  Therefore,  the  convenience  of 
linearity  in  modeling  should  be  abandoned  and  material  and  geometric  nonlinearities  should 
be  accounted  for. 

2.2  Piezoelectric  materials 

Piezoelectricity  is  a  phenomenon  in  which  some  materials  develop  polarization  upon 
application  of  strains^.  This  phenomenon  is  observed  in  materials  that  have  a  non- 
centrosymmetric  crystal  structure.  Examples  of  piezoelectric  materials  are  Rochelle  salt, 
quartz,  and  the  most  popular  one.  Lead  Zirconate  Titanate  or  PZT  (Pb  (Zr,Ti)  O3). 

Piezoelectric  materials  exhibit  a  linear  relationship  between  the  electric  field  and  strains 
for  low  field  values  (up  to  100  V/mm).  However,  the  relationship  is  nonlinear  for  large  fields, 
and  the  material  exhibits  hysteresis^.  Furthermore,  piezoelectric  materials  show  dielectric 
aging  and  hence  lack  reproducibility  of  strains,  ie.,  a  drift  from  zero  state  of  strain  is  observed 
under  cyclic  electric  field  conditions'^. 

2.3  Electrostrictive  materials 

Electrostriction  is  a  second-order  effect  observed  in  all  dielectrics  (centrosymmetric  and 
non-centrosymmetric  crystal  structure)  where  an  applied  uniform  electric  field  induces  a 
strain,  but  a  reversal  of  the  field  does  not  alter  the  strain^.  That  is,  the  strain  is  proportional 
to  the  square  of  the  electric  field.  An  electrostrictive  material  that  is  being  used  increasingly 
is  Lead  Magnesium  Niobate  or  PMN  (Pb  (Mgi/3  Nb2/3)  O3).  These  materials  have  a 
nonlinear  field-strain  relationship,  but  exhibit  very  little,  if  any,  hysteresis®.  Additionally, 
these  materials  have  excellent  zero  strain  reproducibility  in  cyclic  electric  field  conditions®. 
However,  strains  induced  by  electrostriction  are  comparable  to  piezoelectricity  only  at  large 
electric  fields.  But  the  voltage  can  be  reduced  to  the  piezoelectric  level  for  an  electrostrictive 
material  with  a  high  dielectric  constant®.  Table  2.1  summarizes  the  characteristics  of 
piezoelectric  and  electrostrictive  materials. 

2.4  Composites 

The  active  composite  elements  are  available  in  two  forms:  multi-layered  active  materials 
and  multi-phase  active  materials.  The  multi-layered  active  materials  are  stacks  of  active 
materials  wherein  it  is  possible  to  achieve  large  fields  for  moderate  voltage  input  levels®“^^.  A 
typical  multilayer  piezoelectric  actuator  is  shown  in  Figure  2.3.  The  multilayer  active  ceramics 
are  manufactured  using  tape  casting  technology  and  electrodes  are  embedded  in  several 
different  ways®’^®.  The  multi-layer  composite  system  is  also  amenable  to  miniaturization 
as  it  is  possible  to  decrease  the  thickness  of  active  material  layers.  However,  accounting  for 
nonlinearity  in  modeling  is  essential.  Considerable  research  and  development  has  been  carried 
out  in  piezoelectric  multi-layered  ceramics.  Active  vibration  control  using  these  multi-layered 
piezoceramic  configurations  has  been  investigated^®”^’^. 


4 


Table  2.1.  Characteristics  of  Piezoelectric  and  Electrostrictive  Materials. 


Piezoelectric  Materials 

•  At  low  electric  fields,  the  strain-electric  field  relationship  is  linear. 

•  At  large  electric  fields,  the  relationship  is  nonlinear,  and  hysteresis  appears  under  cyclic 
field  conditions. 

•  Exhibits  drift  from  zero  strain  state  under  cyclic  field  conditions  due  to  dielectric  aging. 

•  In  the  linear  regime,  the  piezoelectric  coefiicient  is  constant  and  hence  cannot  be 
electrically  tuned  with  a  bias  field. 

•  Has  been  studied  extensively  in  the  last  few  years,  and  several  analysis  tools  have 
emerged;  experimental  characterization  has  been  done  extensively. 

•  Currently  still  popular  because  of  lower  power  requirements  at  high  frequencies,  relatively 
lower  prices,  and  above  all  an  understanding  of  the  component  integrity  at  low  fields. 

Electrostrictive  Materials 

•  The  strain-electric  field  relationship  is  nonlinear. 

•  Shows  very  little,  if  any,  hysteresis. 

•  Has  excellent  reproducibility  of  zero  strain  state. 

•  Requires  a  high  dielectric  constant  to  generate  a  large  electromechanical  coupling. 

•  Nonlinear  relation  between  the  strain  and  field  can  be  used  to  tune  the  dielectric  constant; 
therefore,  the  electromechanical  coefficient  can  be  tuned  over  a  wide  range,  changing  from 
inactive  to  extremely  active  states. 

•  While  the  constitutive  relationship  has  been  in  existence  in  the  literature  for  quite 
a  while®"®,  modeling  and  analysis  of  the  behavior  of  electrostrictive  materials  under 
various  loading  conditions  have  not  been  carried  out.  Further  understanding  of  reliability 
and  structural  integrity  of  smart  structural  systems  with  electrostrictive  materials  is 
necessary. 


Multi-phase  composite  systems  consist  of  an  active  phase  embedded  in  a  passive  matrix 
(for  example,  PZT  rods  in  a  polymer  matrix^®"^®,  as  shown  in  Figure  2.4).  Different 
connectivity  schemes  are  possible  in  a  two-phase  system.  The  classification  based  on 
connectivity  of  active  and  passive  phases  was  developed  by  Newnham,  et  al.^^.  Literature 
abounds  with  papers  on  piezoceramic  rods  (active)  embedded  in  pol3aner  matrix  (passive), 
forming  a  1-3  connectivity  composite  system^®’^®.  Finite  element  analysis  of  1-3  composite 
systems  has  also  been  carried  out^^’^®.  These  multi-phase  composites  are  most  often  useful 
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when  there  are  conflicting  requirements  on  the  active  material  system.  This  type  of  active- 
passive  phase  combinations  helps  in  achieving  multi-functionality.  However,  for  reliability, 
mechanics  studies  relating  to  active/passive  material  bonding  interactions  are  necessary^^. 
Certain  receptor  cells  in  the  cochlea  of  the  inner  ear  undergo  displacements  flve  orders  of 
magnitude  more  than  what  the  best  piezoelectric  ceramic  can  produce  at  the  same  level  of 
applied  voltage  and  dimensionality^^. 


Figure  2.3;  A  multi-layer  actuator. 


Figure  2.4;  A  schematic  of  a  1-3  connectivity  composite. 


2.5  Issues  in  modeling  and  analysis 

The  mechanics  of  smart  material  systems  involves  coupling  between  electric,  magnetic, 
thermal,  and  mechanical  effects.  In  addition  to  this  coupling,  it  may  be  necessary  to  account 
for  geometric  and  material  nonlinearities.  For  example,  an  electromechanical  transducer  is 
characterized  by  five  important  properties^:  the  resonant  frequency,  acoustic  impedance. 
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mechanical  damping  coefficient,  electromechanical  coupling  coefficient,  and  the  electric 
impedance.  If  nonlinear  electroelastic  equations  are  included  in  the  model,  some  or  all  of  these 
properties  can  be  tuned;  for  instance,  in  an  electrostrictive  material,  the  electromechanical 
coupling  coefficient  can  be  tuned  with  a  bias  field^.  In  order  to  tune  the  first  fundamental 
resonant  frequency  of  the  transducer,  thin  rubber  layers  are  introduced  in  a  multilayer  PZT 
laminate^®.  The  thin  rubber  layers  necessitate  the  use  of  nonlinear  elastic  relations. 

Toupin^^  was  the  first  one  to  consider  the  nonlinearity  in  electroelastic  formulations. 
Knops^®  presented  a  two-dimensional  theory  of  electrostriction  and  solved  a  simplified 
boundary  value  problem  using  complex  potentials.  However,  rotationally  invariant  nonlinear 
thermoelectroelastic  equations  were  derived  by  Tiersten  and  his  coworkers®’^’®,  and  by 
Nelson®.  Tiersten®  has  stressed  the  importance  of  including  nonlinear  terms  in  the 
constitutive  relations,  particularly  at  large  fields.  Joshi®®  has  also  presented  nonlinear 
constitutive  relations  for  piezoceramic  materials. 

To  our  knowledge,  a  fully  coupled,  nonlinear  electrothermomechanical  analysis  has 
not  been  carried  out  so  far.  Analytical  studies  have  largely  concentrated  on  fracture  of 
piezoceramics  and  electrostrictive  ceramics.  Cracks  emanating  from  electrode  junctions, 
debonding  in  active/passive  composite  systems,  and  delamination  in  multi-layered  stacking 
are  some  of  the  analytical  problems  that  have  been  considered  in  the  literature®®’®^ .  Coupled 
analyses  were  carried  out  for  linear  piezoelectric  materials®®. 

During  the  present  research,  a  general  formulation  for  laminated  composite  plates  with 
piezoelectric/magnetostrictive  actuators  and/or  sensors  has  been  developed.  Formulations  of 
both  classical  and  shear  deformation  (first-order  and  third-order)  theories  were  developed,  and 
they  take  into  accoimt  the  thermo-electro-mechanical  coupling,  von  Karman  type  geometric 
nonlinearity,  and  time  dependency.  The  Navier  solutions  for  the  linear  cases  were  also  derived. 
Assumed  kinematics  of  deformation  through  thickness  of  the  plate  in  various  plate  theories 
is  graphically  depicted  in  Figure  2.5. 


Figure  2.5:  Kinematics  of  deformation  of  a  plate  edge  in  various  plate  theories. 
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3.  Theoretical  formulation  using  CLPT 
3.1  Displacements  and  stredns 

The  classical  laminated  plate  theory  (CLPT)  is  based  on  the  KirchhofF  assumption  (see 
Reddy^^-®^)  that  straight  lines  perpendicular  to  the  midsurface  {i.e.,  transverse  normals) 
before  deformation  remain  straight  after  deformation,  they  are  inextensible,  and  rotate  such 
that  they  remain  perpendicular  to  the  midsurface  after  deformation.  These  assumptions 
imply  that  the  transverse  normal  strain  Czz  and  transverse  shear  strains  and  Syz  are  zero. 

Consider  a  plate  of  total  thickness  h  composed  of  N  orthotropic  layers  with  the  principal 
material  coordinates  of  the  kth  lamina  oriented  at  an  angle  Ok  to  the  laminate  coordinate,  x. 
The  xy— plane  is  taken  to  be  the  undeformed  midplane  flo  of  the  laminate,  and  the  z-axis 
is  taken  positive  upward  from  the  midplane.  The  fcth  layer  is  located  between  the  points 
z  =  Zk  and  z  =  Zk+i  in  the  thickness  direction.  Some  of  the  layers  have  the  sole  purpose  of 
actuating  or  sensing  deformation  of  the  laminate  (see  Figure  3.1). 


Figure  3.1;  A  schematic  of  a  laminated  plate  with  imbedded  actuating/sensing  layers. 


In  formulating  the  theory,  we  assume  that  the  layers  are  perfectly  bonded  together. 
Further,  restrict  the  formulation  to  linear  elastic  material  bahavior,  small  strains  and 
displacements,  and  to  the  case  in  which  the  temperature  and  electric  fields  are  given. 

The  Kirchhoff  hypothesis  leads  to  the  displacement  field 

u{x,  y,  z,  t)  =  UQ{x,y,t)  -  z^^ 
v{x,  y,  z,  t)  =  vo{x,  y,t)-  z^^ 

w{x,  y,  z,  t)  =  wo{x,  y,  t)  (3.1) 

where  {uo,vo,wq)  are  the  displacements  along  the  coordinate  lines  of  a  material  point  on  the 
xj/— plane.  Once  the  midplane  displacements  {uq,vo,wo)  are  known,  the  displacements  of  any 
arbitrary  point  {x,y,z)  in  the  3-D  continuum  can  be  determined  using  Eq.  (3.1). 
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The  nonzero  von  Karman  strains  associated  with  the  displacement  field  in  Eq.  (3.1)  are 
gi^enby  ___ 


(3.2a) 
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where  are  the  membrane  strains  and  {sx'x  are  the  flexural  (bending) 

strains.  The  transverse  strains  (sxz^^yz^^zz)  are  identically  zero  in  the  classical  plate  theory. 
Note  from  Eq.  (3.2b)  that  all  strain  components  vary  linearly  through  the  laminate  thickness, 
and  they  are  independent  of  the  material  variations  through  the  laminate  thickness. 

3.2  Lamina  constitutive  relations 

The  linear  constitutive  relations  for  the  A:th  orthotropic  (piezoelectric)  lamina  in  the 
principal  material  coordinates  of  a  lamina  are 

\Qn  Qi2  0  rei-ai  AT]  TO  0  esi  1  f  5i 

<  Cr2  ?  =  Qi2  Q22  0  <62-0^2  at  >  “  0  0  632  {  ^2  ?  (3.3) 

[  (J6  J  ^0  0  Qee.  I  ^6  J  .0  0  0  _  l^sj 

where  are  the  plane  stress-reduced  stiffnesses  and  are  the  piezoelectric  moduli  of 
the  fcth  lamina,  {(Ti^ei^Si)  are  the  stress,  strain,  and  electric  field  components,  respectively, 

referred  to  the  material  coordinate  system  (ori,  0:27  ^3)7  and  0:2  are  the  coefficients  of 

thermal  expansion  along  the  xi  and  X2  directions,  respectively,  and  AT  is  the  temperature 
increment  from  a  reference  state,  AT  =  T  —  Tq.  For  layers  other  than  piezoelectric  layers, 
the  part  containing  the  piezoelectric  moduli  should  be  omitted.  The  coefficients  are 
known  in  terms  of  the  engineering  constants  of  the  fcth  layer: 

gii  =  - ,  Qi2  =  Q22  =  - ,  Q66  =  Gi2  (3.4a) 

and  the  piezoelectric  stiffnesses  are  known  in  terms  of  the  dielectric  constants  and  elastic 
stiffnesses  as 

'0  0  ro  0  dsiY’^UQn  Q12 

0  0  632  =00  ^32  Qi2  Q22  0  (3.46) 

.OOOj  LOOOJ  [0  0  ^66. 

Since  the  laminate  is  made  of  several  orthotropic  layers,  with  their  material  axes  oriented 
arbitrarily  with  respect  to  the  laminate  coordinates,  the  constitutive  equations  of  each  layer 
must  be  transformed  to  the  laminate  coordinates  (x,y,z).  The  transformed  stress-strain 
relations  relate  the  stresses  {axx,(^yy,(^xy)  to  the  strains  {exx,Syy,lxy)  and  components  of  the 
electric  displacement  vector  {SxiSy^E^  in  the  laminate  coordinates 
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where 

Oil  =  Oil  cos*^  9  +  2{Qi2  4-  2(566)  sin^  9  cos^  9  +  Q22  sin"^  9 
O12  =  (Oil  +  O22  —  4066)  sin^  9  cos^  9  +  (^laCsin"^  9  +  cos'^  9) 

O22  =  Oil  sin^  9  +  2{Qi2  +  2(566)  sin^  9  cos^  9  +  Q22  cos^  9 

O16  =  (On  -  O12  -  2066)  sin  9  cos®  9  +  ((5i2  -  O22  +  2066)  sin®  9  cos  9 

O26  =  (On  ~  O12  —  2(566)  sin®  0cos^  +  (O12  ~  O22  +  2066)  sin  ^  cos®  9 

066  =  (On  +  O22  -  2O12  -  2066)  sin®  9  cos®  9  +  (566(sin‘*  9  +  cos^  9)  (3.6) 

and  ai,  a2,  and  fig  are  the  transformed  thermal  coefficients  of  expansion 

q:i  =  «!  cos®  9 -\-a2  sin®  9 
di2  =  sin®  0  +  02  cos®  9 

og  =  2(q;i  —  02)  sin  9  cos  9  (3.7) 

and  Cij  are  the  transformed  piezoelectric  moduli 

631  =  631  cos®  9  +  632  sin®  9  =  Qee  (dsi  cos®  9  +  6.32  sin®  0^ 

632  =  631  sin®  0  +  632  cos®  0  =  (566  (^31  sin®  0  +  d32  cos®  0^ 

636  =  (631  -  632)  sin  0  cos  0  =  Qm (dsi  -  ^32)  sin  0  cos  0  (3.8) 

Here  0  is  the  angle  measured  counterclockwise  from  the  x— coordinate  to  the  xi —coordinate. 
We  define 


where  summation  is  implied  on  repeated  subscript  j  over  the  range  j  =  1, 2,6. 


3.3  Equations  of  motion 

The  equations  of  motion  can  be  derived  using  the  principle  of  virtual  displacements. 
In  the  derivations,  we  account  for  thermal  (and  hence,  moisture)  and  piezoelectric  effects 
only  with  the  understanding  that  the  material  properties  are  independent  of  temperature 
and  electric  fields,  and  that  the  temperature  T  and  electric  displacement  vector  S  are  known 
functions  of  position.  Thus  temperature  and  electric  fields  enter  the  formulation  only  through 
constitutive  equations. 


The  equations  of  motion  of  the  classical  plate  theory  are 
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where  q  is  the  distributed  transverse  mechanical  load,  and  (Nxx^Nyy^Nxy)  the  total  in- plane 
force  resultants  and  (MxxyMyy^Mxy)  the  total  moment  resultants  defined  by 


(7oj  ^1,-^2)  ^ire  the  metss  moments  of  inertia 


(3.13) 
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and 
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Note  that  the  force  and  moment  resultants  are  related  to  the  mechanical,  thermal,  and 
piezoelectric  strains  through  the  constitutive  equations.  We  will  consider  the  laminate 
constitutive  equations  in  the  sequel. 


The  generalized  displacements  and  forces  of  CLPT  are 
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where  Qn  denotes  the  transverse  shear  force  resultant  on  an  edge  with  normal  n  =  (ux^ny). 
In  the  dynamic  case,  it  is  given  by 
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The  force  and  moment  resultants  on  an  edge  with  normal  n  to  those  on  edges  parallel  to  the 
X  and  y  coordinates  by 
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(3.186) 


3.4  Laminate  constitutive  equations 

Here  we  relate  the  force  and  moment  resultants  in  Eq.  (3.13)  to  the  strains.  By  definition, 
we  have 
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where  {N^}  and  {N^}  are  thermal  and  electric  force  resultants 
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Na  is  the  number  of  actuating  layers,  and 
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Similarly,  we  have 
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where  {M^}  and  {M^}  are  thermal  and  electric  moment  resultants 
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Here  Aij  denote  the  extensional  stiffnesses,  Dij  the  bending  stiffnesses,  and  Bij  the  hending- 
extensional  coupling  stiffnesses 
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Dij  =  \'EQf(4+j-4) 


(3,26) 


Note  that  Q’s,  and  therefore  yl’s,  B’s,  and  £)’s,  are,  in  general,  functions  of  position  {x,y). 
Equations  (3.19)  and  (3.23)  can  be  written  in  a  compact  form  as 
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Suppose  that  the  temperatvire  and  electric  fields  vary  linearly  within  feth  layer 

AT  =  T^{x,y,t)i;!l{z)  +  r'=(x,y,  t)rpl{z)  (3.28a) 

Sz  =  S^{x,y,t)'il;'l{z)  +  S^{x,y,t)i)^{z)  (3.286) 

where  are  the  linear  interpolation  functions  of  the  feth  layer 
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and  Ti  and  T2  denote  the  temperature  at  the  bottom  and  top  surfaces  of  the  fcth  layer. 
Similar  notation  is  used  for  the  electric  field.  Then  the  thermal  and  electrical  forces  and 
moments  can  be  evaluated  as 
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On 

012 

O16 

(fc)  I 

rd}i 

012 

022 

O26 

< 

I  ^32 

.016 

026 

Oee. 

I  ^36 

(k) 


Al 


{T^^\  +  T^^l;l)dz={  ^  (3.30) 

Al 


(k) 


Bl 


(TiVi  +  T^i^2>dz  =  {  Bj  K3-31) 


(fc) 


{e^tl^^,+£^xP^)dz  = 


kj.k\J _ 


(fc) 


Bi 


(£!^tpt  +  £2^2)^dz={  B^  k3.33) 


Bi 


^  fc=lj=l,2,6 

Bl  =  lE  E  [It  {hk  +  3zfc)  +  T2  {2hk  +  ^Zk)^  hk 


A:=:li=l,2,6 
are  the  thermal  stiffnesses,  and 

Na 

2 


1  i»a 

^  fc=l  i=l,2,6 

=  E  Q\f^3^[sHhk  +  ^Zk)  +  £H‘2hk  +  Szk)]hk 


(3.34a) 

(3.346) 

(3.35a) 

(3.356) 


k=l  i=l,2,6 
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are  the  piezoelectric  stiffnesses,  and  Na  <  N  is  the  number  of  actuator  layers.  If  the  electric 
field  intensity  is  constant  across  each  lamina,  then  one  may  approximate  it  as 


CK  _  ck  _ 

—  ^2  — 


(3.36) 


where  Vfe  is  the  applied  voltage  across  the  fcth  layer  and  hk  is  the  thickness  of  the  layer. 

The  stress  resultants  (iV’s  and  M’s)  are  related  to  the  displacement  gradients, 
temperature  increment,  and  electric  field  by  the  relations 


^11  ^12  Ai6 
Ai2  A22  A26 
Aie  A2e  Aqq 


Bn  B\2  Biq 

B\2  B22  ^26 

Bie  B26  Bee 


dug  \  X  ( 
dx  '  2  0a;  j 

dvg  \  X  ( ^^^0 
dy  ^  2  \  dy  ) 
dug  I  dvg  _j_  dwg  dwg 
dy  dx  dx  dy 


dug  \  X  ( 
dx  '  2  \  dx  J 

dvg  i  X  { 

dy  2  \  dy  J 

dug  I  dvg  I  dwg  dwg 

dy  '  dx  *  dx  dy 


Bn  Bi2  Bie 
Bi2  B22  B2e 
Bie  B2e  Bee 


Dn  Di2  Die 
Di2  D22  D2e 
Die  D2e  Dee 


(3.37) 


(3.38) 


4.  Theoretical  formulation  of  TSDT 


4.1.  General  comments 


The  classical  laminate  plate  theory  and  the  first-order  shear  deformation  theory  are  the 
simplest  equivalent  single-layer  theories,  and  they  adequately  describe  the  kinematic  behavior 
of  most  laminates^^.  The  third-order  theory  of  Reddy®^’^^’^  represents  the  plate  kinematics 
better,  does  not  require  shear  correction  factors,  and  yields  more  accurate  interlaminar  stress 
distributions.  Hence,  the  generalized  third-order  shear  deformation  laminate  theory  (TSDT) 
of  Reddy^^,  which  contains  the  first-order  shear  deformation  theory  (FSDT)  as  a  special  case, 
is  used  to  develop  the  governing  equations  of  laminated  plates  with  actuating  and/or  sensing 
layers. 


4.2.  Displacement  field  and  strains 

The  Reddy  third-order  plate  theory  is  based  on  the  displacement  field 

u{x,y,z,t)  =  uo{x,y,t)  +  z<}>x{x,y,t)  -  ciz^ 

v{x,  y,  z,  t)  =  vo{x,  y,  t)  -h  z<f>y{x,  y,  t)  -  ciz^  (^(t>y  -H 
w{x,y,z,t)  =wo{x,y,t) 


(4.1) 
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where  (uq,vo,wo)  and  {(f>x,<f>y)  have  the  same  physical  meaning  as  in  the  first-order  theory; 
they  denote  the  displacements  and  rotations  of  transverse  normals  on  the  plane  z  =  0, 
respectively.  Then  the  displacement  field  of  FSDT  is  obtained  by  setting  ci  =  0,  and  for  the 
Reddy  third-order  theory  we  set  ci  =  4/(3h^). 

Substitution  of  the  displacements  (4.1)  into  the  von  Karman  nonlinear  strain- 
displacement  relations  yields  the  strains^^ 


where 


(4.4) 

(4.5) 

(4.6) 


4.3.  Constitutive  relations 

In  addition  to  the  constitutive  relations  in  Eq.  (3.5),  we  have  the  following  requations 
for  the  transverse  shear  stresses 

/  \  1  /  \  4-  [  ®  \  fA  7\ 

{cra^z}  [Q45  Q55J  Wzr  [^15  ^25  Oj  |  g  J 

where 

O44  =  Qu  cos^  6  +  Q55  sin^  6 

O45  =  (Q55  “  Q44)  cos  ^  sin  0 

O55  =  Q55  cos^  6  +  Qu  sin^  0  (4.8) 

ei4  =  (ei5  —  624)  sin  6  cos  9 
624  =  C24  cos^  6  +  ei5  sin^  6 
615  =  ci5  cos^  9  +  624  sin^  9 
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Also,  we  have 


625  =  (ei5  —  624)  sin  0  cos  0 


611  =  eii  cos^  9  +  €22  sin^  9 
622  =  611  sin^  9  +  €22  cos^  9 

612  =  (fn  -  622)  sin  ^  cos  0 


614  615 

624  625 

0  0 


di4  di5 
^24  d25 

0  0 


Q44  Q45 

Qa5  Q55 


(4.9) 


(4.10) 


(4.11) 


4.4.  Equations  of  motion 

The  equations  of  motion  may  be  developed  using  the  dynamic  version  of  the  principle  of 
virtual  displacements. 


where 


dNxj.  .  dN. 


+ 


xy  _ 


dx 

dNxy  .  dNi 


dy 


Iquq  +  J\^x  —  C1/3 


dx 

dQx 


+ 


yy 


% 

dQy 


=  IqVQ  +  J\^y  —  Cl  J3 


dwp 

dx 

dwo 


dy 
dwo 


d  ,  dwo 

~E - 1"  - h  Nxy-^ 

OX  ay  ox  ox  ay 


)  +  l.(N  ^  +  N  — ° 
’^dy^^^^  dx  dy 


d^Px: 


^  ^dilo  ,  dvo  d^x  ,  d<j)y^ 


+  Ci( 

+  Ci 

dx 

dMxy  .  dM^ 
dx 


+  2 


d^^Px^  d'^P 


+  ■ 


yy 


dx^  dxdy  dy^ 


)  +  q  =  IqWo  -  c\h{ 


d^wo  di^WQ 
dx?'  dy’^ 


+  • 


dM, 


xy 


dy 


+ 


yy 


dy 


—  Qx  =  JiUq  +  K2<I>x  —  Cl  J4 


Qy  —  d\VQ  ”i”  K^^y  Cl .74 


dwp 

dx 

d'WQ 

dy 


Ca  Sc^Roi 

p^''\zy  dz  {i  =  0,1, 2,  ■■■,6) 

I— -1  Jzu 


—  ^a(3  3!)9?  Qa  —  Qc 

N 


Ji=  li  —  Ci/i4-2,  K2  =  h  —  2ci/4  +  cf/e,  C\  = 


and  {Pxx^PyyjPxy)  ^nd  (RxyRy)  denote  the  higher-order  stress  resultants 


yy 

o 

xy 


(s)-C(:n} 


z^dz 


(4.12) 

(4.13) 


(4.14) 

(4.15) 

(4.16) 

(4.17) 
(4.18a) 

(4.185) 

(4.19) 
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The  primary  and  secondary  variables  of  the  theory  are 

Primary  Variables  :  Un,  Ug,  wq,  <i>n,  <t>s 

on 

Secondary  Variables  :  Nnn,  Nns,  Pnn,  Mnn,  Mns 
where  has  the  form  similar  to  Mnn  [see  Eqs.  (2.18o,6)]  and 


hn  — 


Cl 


/  dPxx  .  dPxy  \  fdPxy  dPyy\ 

V  9a:  dy  )  ^  V  9a:  dy  ) 

"h  Ji^x  CiIq  ^  nx  ~t"  ^-^3^0  Cl Iq  ^ 


+ 


dPn 


(^Qx'^x  “h  Qy'C^y}  ”i”  ^xxy  ^xy}  T  Ci 


V={Nxx^+Nxy 


dwQ 


+  Ny 


9n;o\ 


xy  -r-j/y  Qy  jny 


(4.20a) 

(4.206) 


(4.21) 

(4.22) 


The  stress  resultants  are  related  to  the  strains  by  the  relations 


[B]  [E] 

[D]  [F] 
[F]  [F]J 


(Aij,Bij,Dij,Eij,Fij,Hij)  =J2  f  Qif  {l,z,z^,z^,z'^,z^)  dz 

1 _ 1  J  Zk 


fc==:l 
N 


{Aij,Dij,Fij)  =  Y^  {lyz\z^) 

k=l  ''^»= 


dz 


(4.23) 


(4.24) 


(4.25a) 


(4.256) 


The  stiffnesses  Aij,  Dij  and  Fij  are  defined  for  i,j  =  1,2,6  as  well  as  i,jf  =  4,5.  The 
stiffnesses  Bij^  Eij  and  Hij  are  defined  only  for  i^j  =  1,2,6.  The  coefficients  Aij^Bij^  and 

Dij  were  given  in  terms  of  the  layer  stiffnesses  and  layer  coordinates  Zk-^i  and  in  Eq. 
(2.26).  The  higher-order  stiffness  coefficients  are  defined  by 

^  fc=l 

Pij  =  ^  [i.^k+\f  -  {Zkf] 

Hij  =  ^  E  Qif  (4-26) 

'  fe=i 

Note  that  the  stiffnesses  Eij,  Fij  and  so  on  of  the  third-order  theory  involve  fourth  or  higher 
powers  of  the  thickness,  and,  therefore,  they  are  expected  to  contribute  httle  to  thin  laminate 
solutions.  Even  for  moderately  thick  laminates  the  contribution  can  be  small. 
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5.  Sensor  sind  actuator  equations 

5.1.  Electro- mechanical  coupling 

The  electro-mechanical  coupling  is  a  two-way  coupling.  The  effect  of  applied  electrical 
field  on  stresses  is,  as  indicated  by  Eq.  (2.5),  is  analogous  to  the  effect  of  the  temperature 
field.  The  components  Dj  of  electrical  displacement  vector  are  related  to  the  components  of 
strains  and  electrical  field  by 

rO  0  Ten  0  0  I  f  5i 'I 

i  D2  ^  =  0  0  0  {  £2  >  +  0  e22  0  I  £2}  (5.1) 

[  T>3  J  .631  632  0.  leej  [o  0  ^33]  I  ^3  J 

where  Sij  are  piezoelectric  stiffnesses  and  €ij  are  dielectric  coeflScients.  The  transformed 
equations  are 

ro  0  o]('')r£xx'i  \hi  ei2  0  r  4 'I 

<  Dy  >  =  0  0  0  I  eyy  >  -  €12  €22  0  I  ^2/  7  (5.2a) 

\Dz]  [e3i  632  633]  l7xyJ  LO  0  €33]  \£z] 

Equation  (5.2a)  is  modified  to  account  for  the  transverse  shear  strains  in  TSDT 

£xx  "I  614  ®15  f  -i  ^11  ^12  0  f  1 

Syy  >  +  624  625  1  f  ^12  ^22  0  S  /  (5-25) 

7xJ  Lo  0]’-^='^'’  [0  0  633]  UJ 

This  electromechanical  equation  provides  the  starting  point  for  the  derivation  of  the 
sensor  and  actuator  equations^^’^^’^®. 

5.2.  Sensor  equations 

According  to  the  Gauss  law,  the  closed  circuit  charge  measured  through  the  electrodes 
of  the  A;th  layer  is 

Qk{t)  =  l  f  Dxdxdy  -f  f  D.dxdy  (5.3) 

[JSj{z=Zh)  JSi(z=Zk+i)  J 

where  Si  ^  SbO  St  is  the  intersection  of  the  electrode  surfaces  on  both  sides  of  the  lamina 
(see  Figure  5.1)  and  subscript  ‘s’  denotes  sensor.  The  electrodes  are  assumed  to  be  placed 
on  the  transverse  surfaces  with  the  poling  direction  z.  Hence,  only  the  component  Dz  of 
the  electric  displacement  vector  is  nonzero.  When  Eq.  (5.2)  is  applied  to  sensors  where  the 
converse  piezoelectric  effect  is  negligible  and  the  external  applied  charge  is  zero,  the  sensed 
charge  can  be  calculated  from  Eq.  (5.2)  as 

{^xx  I 

Syy  >  (5.4) 

Ifxy  J 

The  total  charge  in  the  laminate  is  calculated  by  summing  over  the  number  of  sensor  layers, 

Qs{-t)  =  lf]  I  D,dxdy+  f  D,dxdy  (5.5) 

2  \  JSi{z=Zk)  JSi{z=Zk+i) 


Da,]  r  0  0  0  1  r 

T>j,  I  =  0  0  0  < 

Dz  )  L®31  632  633]  [ 


Substituting  Eq.  (5.4)  into  (5.5),  we  obtain 


N, 

k=l 

where  =  0.5(2^  +  Zk+i) 


Qsit)  =  ^  E  [e3i4i  +  e32e5,?  +  eseT^^  +  4  (e3i4x  +  6324?/  +  63674'’)] 


(5.6) 


Figure  5.1:  Laminate  with  surface  electrodes  {St  ^surface  electrode  in  the  top  surface, 
Sb  =  surface  electrode  in  the  bottom  surface. 


We  note  that  Eq.  (5.6)  was  obtained  by  setting  the  externally  applied  fields  and 
to  zero,  because  when  the  surface  electrodes  are  short-circuited  on  both  sides  of  the 

piezoelectric  lamina,  the  terms  containing  the  externally  applied  electric  field  intensities 
(2) 

and  Sz  ^  should  be  set  to  zero.  The  resulting  equation  is  the  closed-circuit  charge  sensor 
equation^  which  relates  the  output  signal  to  the  plate  deformation. 

The  current  I{t)  on  the  surface  of  the  sensor  is  given  by 

(5.7) 

When  the  sensors  are  used  as  strain  rate  tensors,  the  current  can  be  converted  to  the  open 
circuit  sensor  voltage  ouput  Vs 


Vs{t)  =  Gclit)  = 


(5.8) 


where  Gc  is  the  gain  of  the  current  amplifier. 

The  sensor  forces  {N^}  and  moments  {M^}  (and  {P^})  are  then  evaluated  by  using 
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in  Eqs.  (3.32)  and  (3.33).  Here  hs  denotes  the  thickness  of  the  sensor  layer.  The  forces  and 
moments  will  be  determined  in  terms  of  of  the  strain  rates  in  the  sensor  layers. 


5.3.  Actuator  equations 

The  actuator  equations  for  piezoelectric  actuator  can  be  derived  using  induced  strain 
actuation.  We  assume  that  the  piezoelectric  actuator  layer  does  not  have  any  applied  stress. 
Then  Eq.  (5.1)  gives  the  strains  developed  by  the  electric  field  on  the  actuator  layer 

{£}.  =  iQi-‘ief  {«}  =  ie|-‘iai’'i<i|’'{f}  (5.10) 

or 

=  (5.11) 

The  stresses  due  to  the  actuator  strains  are 

W}a^[Q]{e}a^[Qm'^{£}  (5.12a) 

or 

<  -  Qijdkjek  {ij  =  1, 2, 6;  =  1, 2, 3)  (5.12ft) 

The  stress  vector  can  be  used  to  generate  the  actuator  forces  {N^}  and  moments  {M^}. 

If  the  voltage  applied  to  the  actuator  in  the  thickness  direction  (only)  is  V^,  then 


where  ha  is  the  thickness  of  the  actuator  layer.  Then  we  have 

=  QijdsjVa 

Mf  =  QiidzjVazt 


(5.13) 


(5.14a) 

(5.146) 


where  2:“  is  the  distance  from  the  midplane  to  the  center  of  the  actuator  layer. 


6.  The  Navier  solutions  of  CLPT 
6.1.  Preliminary  comments 

The  equations  of  motion  (3.10)-(3.12)  can  be  expressed  in  terms  of  displacements 
(uo,vo,wo),  temperature  field  (T^,T^),  and  electrical  field  by  substituting  for  the 

force  and  moment  resultants  from  Eqs.(3.37)  and  (3.38).  For  homogeneous  laminates  (i.e., 
for  laminates  with  constant  j4’s,  B's,  and  D’s)  and  under  the  assumption  of  small  strains, 
displacements  and  rotations,  the  equations  of  motion  (3.10)-(3.12)  can  be  expressed  as®^ 


cii  ci2  cisll'uo'l  r^nii  0  mia'irfio'j  j'0'|  f /f 'j  f /f 'j 

C12  C22  C23  <  >  +  0  77122  "»23  <  %  ?  =  S  0  >  +  <  /J  >  +  <  /f  >  (6.1) 

Cl3  C23  C33J  [woJ  L»^13  ^23  ^33  J  (  Wq  J  UJ  I  fs  )  I J 
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where  the  coefficients  Cij  and  rriij  are 


Cii  —  Aiid^  +  2AiQdxdy  +  A^^dy 

Ci2  =  Aisd^  +  (Ai2  +  AQQ)dxdy  +  A2sdy 

ci3  =  —  +  SBied^dy  +  {B12  +  2Be6)dxdy  + 

C22  =  Aeed^  +  2A2&dxdy  +  A22dy 

C23  =  —  [-Biedj  +  {Bi2  +  2Be6)d^dy  +  ZB2&dxd^  +  -822^^ 

C33  =  -Diic^  +  ADied^dy  +  2{D\2  +  2D66)d‘^dy  +  4D2edxdy  +  D22d!^ 

mn  =  -7o4  .  ”^13  =  hdxdt  ,  m22  =  -lodf 

wi23  =  hdyd^  ,  m33  =  /odf  -  hdtidl  +  d^) 


d^xid^yi  d\  denote  the  differential  operators 


cf 

®  dx^ 


(f  =  —  rf  =  — 
y  dy^  ’  ‘  dP 


=  1,2,3,4) 


and  ff  and  are  equivalent  generalized  thermal  and  piezoelectric  forces 


/r  = 


dNix  .  dNjy 


dx 


+  ■ 


,  fi  = 


dN^  dN^ 


fT  _  _  [  ^Ix 

-  y  dx^ 


dy  ^  dx 


+ 


yy 


dy 


XX  ^ 2'' 


+ 


f[  = 


dNP  dNP 


dx 


+ 


xy 


dy 


dxdy  dy^ 


dx 


+ 


dy 


dy^  ) 


f[  =  - 


_ XX  I  fy 

dx^  dxdy 


(6.2) 


(6.3) 


(6.4) 


(6.5) 


In  this  section,  we  develop  the  Navier  solutions  of  Eq.  (6.1)  for  simply  supported 
rectangular  laminates  with  actuating  and/or  sensing  layers.  In  the  Navier  method  the 
generalized  displacements  and  applied  loads  are  expanded  in  a  double  trigonometric  series 
in  terms  of  unknown  parameters.  The  choice  of  the  functions  in  the  series  is  restricted  to 
those  which  satisfy  the  boundary  conditions  of  the  problem.  Substitution  of  the  displacement 
expansions  into  the  governing  equations  should  result  in  a  unique,  invertible,  set  of  algebraic 
equations  among  the  parameters  of  the  expansion.  Otherwise,  the  Navier  solution  cannot  be 
developed  for  the  problem. 

The  Navier  solutions  can  be  developed  only  for  two  classes  of  laminates.  The  first  class 
has  the  following  stiffness  characteristics: 

Aiq  =  ^26  =  Biq  =  ^26  =  £^16  =  ^26  =  0 

Al  =  Bj  =  A^  =  Bi  =  0  (6.6) 

and  the  second  one  has 

Aie  =  j426  =Bii  =  Bi2  =  B22  =  Bee  =  Bie  =  D26  =  0 

A^  =  Bj  =  A^  =  Bi  =  0,  Ji=0  (6.7) 
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These  restrictions  necessarily  require  that  the  actuating  and  sensing  layers  should  be  placed 
such  that  the  laminate  stiffnesses  satisfy  the  above  requirements.  The  stiffness  characteristics 
in  Eq.  (6.6)  are  satisfied  by  antisymmetric  cross-ply  laminates,  while  those  in  Eq.  (6.7)  by 
antisymmetric  angle-ply  laminates^^. 


6.2.  Antisymmetric  cross-ply  plates 

The  simply  supported  (SS-1)  boundary  conditions  used  for  antisymmetric  cross-ply 
laminates  are 


«o(aJ,0,f)  =  0,  ^^{x,0,t)  =  0,  uo{x,b,t)  =  0,  ^^{x,b,t)  =  0 


dx 


M0,y,t)  =  0,  ^(0,y,t)  =  0,  vo{a,y,t)^0,  ^{a,y,t)  =  0 

wo{x,0,t)  =  0,  ■wo{x,b,t)  =0,  wo{0,y,t)  =  0,  WQ{a,y,t)=0 
Nxx{0,y,t)  =  0,  Nxx{0',y,t)  =  0,  Nyy{x,0,t)  =  0,  Nyy{x,b,t)  =  0 
Mxx{0,y,t)  =0,  Mxx{a,y,t)  =  0,  Myy{x,0,t)  =  0,  Myy{x,b,t)  =  0 
where  a  and  b  are  the  planform  dimensions  of  the  plate. 


(6.8) 


The  displacement  boundary  conditions  in  (6.8)  are  satisfied  by  assuming  the  following 
form  of  the  displacements: 


oo  oo 

uo(x,y,t)  =  ^  ^  Umnit)  cos  ax  sin  Py  (6.9a) 

n=l  m=l 
oo  oo 

VQ{x,y,t)  =  EE  VTnn{t)sinax  cos  Py  (6.96) 

n=l  m—1 
oo  oo 

wo{x,y,t)  =  E  Wmn{t)^nax  sin  Py  (6.9c) 

n=l  m=l 


where  a  =  mn/a  and  p  =  mr/b  and  {Umn,Vmn,Wmn)  are  coefficients  to  be  determined. 
Similarly,  the  applied  loads  are  expanded  in  double  sine  series  as 

OO  OO 

q{x,y,t)  =  E  E  Qmn{t)^nax  sin  Py  (6.10a) 

n=l  m=l 
oo  OO 

=  E  E(^^n}sinax  sin;0y  (6.106) 

n=l  m=l 
oo  oo 

{M^}  =  E  E(-^^n}sino:a;  sin  py  (6.10c) 

m=l 
oo  oo 

=  E  E{^mn}sinax  sin  Py  (6.10d) 

n==l  m=l 

n=l  m=l 


(6.10e) 
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where  the  coefficients  Qmni  for  example,  are  given  by 

4  rb 


4  /*o  ra 

Qmnit)  = —r  /  /  q{x,y,t)  sincca:  sin Py  dxdy 

ao  Jo  Jo 


(6.11) 


Substitution  of  Eqs.  (6.9)-(6.11)  into  Eq.  (6.1)  yields  for  any  m,  n,  x,  and  y,  the  following 
differential  equations  in  time 


'cii 

C12 

Cl3‘ 

C12 

C22 

C23 

.^13 

C23 

C33. 

mil  0 

0  17122 

l-ha  -hP 


■ 

Umn  ) 

< 

Vmn 

1 

[Wran] 

+ 


-Iia 

-hP 

^33 

-mLn 


where  Cij  and  mij  are 

cii  =  {Auo?'  +  AeoP'^),  ci2  =  (^12  +  Aoo)oiP 

ci3  =  —Biio?  —  (jBi2  +  2Bee)aP'^,  622  =  {Aooa^  +  A22P^) 

C23  =  —{Bi2  +  ‘^Boo)oC^P  —  B22P^ 

033  =  +  2(Z)i2  +  2Doq)o?P^  +  D22p^ 

mil  =  -fo,  ^22  =  .fo,  ^33  =  io  +  h{oP  +  P"^) 


(6.12) 


(6.13) 


Equation  (6.12)  must  be  modified  to  include  the  sensor  and  actuator  equations.  If  the 
plate  has  only  actuators,  then  the  modification  is  simple;  the  actuator  forces  {N^}  and 
moments  {M^}  in  the  above  equation  must  be  replaced  using  Eqs.  (6.14a,b).  One  must 
be  careful  to  include  them  only  for  the  actuator  layers,  and  for  other  layers  one  must  set 
{ATf }  =  {M'-}  =  {0}. 

6.3.  Antisymmetric  angle-ply  plates 

The  simply  supported  (SS-2)  boundary  conditions  for  the  antisymmetric  angle-ply 
laminates  are 

«o(0,y,t)  =  0,  uo{a,y,t)  =  0,  vo{x,0,t)  =  0,  vo{x,b,t)  =  0 

— (x,0,t)  =  0,  — (x,6,t)  =  0,  — (0,y,t)  =  0,  —{a,y,t)  =  0 

wo{x,0,t)  =  0,  wo(x,b,t)  =  0,  wo(0,y,t)  =  0,  wo(a,y,t)  =  0 
-^xy(0,y,t)  =  0,  Nxy{a,y,t)  =  0,  Nxy{x,0,t)  =  0,  Nxy{x,b,t)  =  0 
Mxxi0,y,t)  =  0,  Mxx{a,y,t)  =  0,  Myy{x,0,t)  =  0,  Myy{x,b,t)  =  0  (6-14) 

The  displacement  boundary  conditions  in  (6.14)  are  satisfied  by  assuming  the  following  form 
of  the  displacements 


uoix,y,t)  Urnn{t)  sinax  cospy 

n~l  m=:l 
00  00 

vo{x,y,t)  =  EE  Vmn{t)  cosax  Sin Py 

n=l  m—l 
00  00 

wo{x,y,t)  =Y^Y1  sin  ax  sin/?y 

n=l  m=l 


(6.15) 

(6.16) 


(6.17) 
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Substitution  of  Eqs.  (6.15)-(6.17)  and  (6.10a-d)  into  Eqs.  (6.1)  yields 


’Cll 

C12 

Cl3 

(  Umn  \ 

'mil 

0 

0  ■ 

C12 

C22 

C23 

\  Vmn  r  “1“ 

0 

rh22 

0 

C23 

C33  +  S33 . 

{Wmn] 

0 

0 

(  0  ]  (  ]  (  -mLn 

=  S  0  >  +  <  >  +  <  -aN^^ 

\Qmn]  I  J  I 

where 


cii  =  A\\a^  +  Aqq0^ ,  ci2  =  {Ai2  +  A^)aj5 

ci3  =  — (SBieo:^  +  B2q0^)P,  C22  =  AQQa^  +  ^22/?^ 

C23  =  —{Biea^  +  3B26l3^)oi 

C33  =  Diia^  +  2(Di2  +  2Dqq)oP‘0^  +  D22^ 

=  .fO)  ^22  =  Iq,  ^33  =  To  +  l2{c?  +  0^) 


(6.18) 


(6.19) 


7.  The  Navier  solutions  of  TSDT 
7.1.  Preliminzury  comments 

The  five  equations  of  motion,  Eqs.(4.12)-(4.16),  admit  the  Navier  solutions  for  simply 
supported  antisymmetric  cross-ply  and  angle-ply  laminates®^’®®*®^.  For  antisymmetric  cross- 
ply  laminates  the  following  stiffnesses  are  zero: 


^16  =  ^26  =  -^45  =  ^16  =  ^26  =  DyQ  =  D26  =  Ji  =  0 


Ei6  =  £'26  =  £i6  =  £26  =  =  H2&  —  £45  =  £45  =  h  =  h  =  I7  =  0  (7.1) 

and  for  antis3anmetric  angle-ply  laminates  the  following  stiffnesses  are  zero: 


.^16  =  .<426  =  -<445  =  £11  =  £12  =  £22  =  £66  =  £l6  =  ^26  =  II  =  0 


£11  =  £12  =  £22  =  Bee  =  £i6  =  £26  =  Hie  =  £26  =  £45  =  £45  =  H  =  H  =  h  =  0  (7.2) 

The  Navier  solutions  of  the  linear  equations  for  these  two  cases  are  presented  next. 

7.2.  Antisymmetric  cross-ply  plates 

The  SS-1  boundary  conditions  for  the  third-order  shear  deformation  plate  theory  are 
uo(aJ,0,t)  =  0,  (f>x{x,0,t)  =  0,  uo{x,b,t)  =  0,  (f>x{x,b,t)  —  0 

vo(0,y,i)  =  0,  <t>y(0,y,i)  =  0,  vo{a,y,t)  =  0,  <f)y{a,y,t)  =  0 
wo(x,0,t)  =  0,  wo(x,b,t)  =  0,  wo(0,y,t)  =  0,  wo(a,y,t)  =  0 
Hxx(0,y,t)  =  0,  iVxa:(o,i/,  i)  =  0,  Nyy(x,0,t)  =  0,  JVyy(x,b,t)  =  0 

■^xx(0,y,t)  =  0,  Afii(a,j/,t)  =  0,  Myy(x,0,t)  =  0,  Myy(x,b,t)  =  0  (7-3) 
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The  boundary  conditions  in  Eq.  (7.3)  are  satisfied  by  the  following  expansions: 

oo  oo 

uo{x,y,t)  =  XI  51  Umn{t)  cos  ax  sin  Py  (7.4o) 

n=l  m=l 
oo  oo 

^0  (a:,  t/,  <)  =  X  X  (*) 

n=l  m=l 
oo  oo 

wq{x, y,t)  =  X  X  W^mn(i)sinaa;sin/?y  (7.4c) 

n=l  m=l 
oo  oo 

=  L  E  Xmn  (t)  COS  ax  sin  Py  {7  Ad) 

n=l  m=l 
oo  oo 

<py  (x,  ?/» ^)  =  X  X  ^rnn{t)  sin  ax  cos  Py  (7.4e) 

Ti=l  m=l 

The  transverse  load  q,  thermal  forces  and  moments  {{N^},{M^},{P^}),  and  electrical 
forces  and  moments  {{N^},  {M^},  {P^})  are  also  expanded  in  double  Fourier  sine  series 
[see  Eqs.(6.10a-e)]. 

Substitution  of  Eqs.(7.4a-e)  and  applied  load  expansions  into  Eqs.(4.12)-(4.16),  we  obtain 
a  5  X  5  system  of  the  following  differential  equations  in  time 

'111  Si2  Sl3  514  515'  (  Umn  "j  'din  0  0  0  ^  '  f  1 

512  522  523  524  525  0  77122  0  0  0  Vmn  I 

513  523  533  S34  S35  <  Wmn  >  +  0  0  77133  1^34  77I35  <  Wmn  / 

514  524  534  S44  S45  Xmn  0  0  77I43  77I44  0  Xmn  I 

.  515  525  535  S45  S55 .  ^  Yfjin  >  0  0  77153  0  77155  J  V  ^mn  ) 

{  0  4  r«iV^nl 

0  0X2mn 

~  ^  Qmn  ^  ^  ^  ^  ^  (^*^) 

0  ^^Imn 

I  0  J  [pMl^^)  [pMi^^) 

where  Sij  and  ihij  are  defined  by 

5ii  =  .diia^  +  AqqP"^,  si2  =  (>li2  +  AQ^)aP 
si3  =  —Cl  I^jEhq:^  +  (P12  +  2J5'66)^^]  « 

5i4  =  Piio:^  +  BeeP^,  sis  =  (.B12  +  Bm)aP 

522  =  -466  Cl^  +  A22P^,  hi  =  5i5 

523  =  —Cl  [p22/?^  +  (P12  +  2^66)“^]  P,  hs  =  Beeoi^  +  B22P‘^ 

533  =  ^55^^  +  ^44/5^  +  <?x  [HuCi^  +  2{Hu  +  2Hqq)oPp'^  +  H22p^] 

534  =  -A55Q;  —  Cl  |^Ai«^  +  (A2  +  2FeQ)aP^ 

535  =  -444/5  -  Cl  +  (A2  +  2Fse)a'^P 

S44  =  -455  +  -Diio:^  +  DeeP^ ,  S45  =  (P12  +  D^)aP 
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555  =  ^44  +  DqqO^  +  D220^  (7.6a) 

mil  =  h,  ^22  =  h,  ^33  =  7o  +  cf/e  [a^  +  0^^  ,  77134  =  -ci  J4a 
^35  =  -CiJ4/3,  77144  =  ^2,  77155  =7^:2  (7.65) 

Aij  =  Aij  c\Dij^  Bij  =  Bij  c\Eij^  Dij  —  Dij  (i, jf  =  1, 2,6) 

Fij  =  F ij  cxHijy  Aij  =  Aij  C\i)ij  =  2ciZ)ij  +  Cii^ij  (7,  j  =  1, 2, 6) 

Dij  =  Aj  -  Cl  A  =  Dij  -  2ciFij  +  clHij  {iJ  =  1,2, 6) 

3ciZ7^j  =  Aij  Oc^Dij  ^c^Fij  (J'^j  ~  ^j^)  (7.6c) 

The  ordinary  differential  equations  (7.5)  in  time  can  be  solved  for  transient  response  using 
the  Newmark  integration  procedure. 

7.3.  Antisymmetric  angle- ply  plates 

The  SS-2  boundary  conditions  for  the  third-order  shear  deformation  plate  theory  are 
^0(0,  y,  t)  =  0,  uo{a, y,  t)  =  0,  vo{x, 0,  t)  =  0,  vo{Xy  5,  t)  =  0 
4>x{xy0,t)  =  0,  ^x(^,5,0  =  0?  <5^2/(0,y,t)  =  0,  ^y{a,y,t)=0 
wo{x,0,t)=0,  wo{x,b,t)  =  0,  7/;o(0,7/,t)  =  0,  wo{a,y,t)—0 
Nxy{0,y,t)  =  0,  Nxy{a,y,t)  =  0,  Nxy{x,Q,t)  =  0,  Nxy{x,b,t)  =  0 
'^17a:x(07  y?  ^)  “  0>  ■^l^xx(c5'j  yj  ^)  —  0,  A7y2/(ir,  0,  t)  =  0,  A7^y(x,  6,  t)  =  0  (^•^) 

The  simply  supported  (SS-2)  boundary  conditions  in  Eq.(7.7)  are  satisfied  by  the 
displacement  expansions 

oo  oo 

uo{x,y,t)  =  x:  Umn{t)  cos  ax  sin  Py  (7.8a) 

n=l  771=1 

oo  oo 

vo{x,y,t)  =  EE  Vmn{t)  sin  ax  cos  py  (7.86) 

77=1  771=1 

and  the  expansions  of  (wq,  <i>x,  4>y)  in  Eqs.(7.4c,d,e).  Substituting  these  expansions  along  with 
the  load  expansions  into  Eqs.  (4.12)-(4.16),  we  obtain  equations  of  the  form  in  (7.5) 

(S){A}  +  1M|{A}  =  {f }  (7.9) 

with  the  following  definition  of  the  coefficients; 

sii  =  +  AeeP^ ,  su  =  {A\2  +  Aee)aP,  I13  =  — ci  (SEiea"^  +  ^26/3^)  P 

si4  =  2BieaP,  sjs  =  Bi^a^  +  B2qP^,  S22  =  j4660:^  +  A22P‘^ 

§23  =  —Cl  ^Eieo;^  +  3^26/3^)  a,  §24  =  §15,  §25  =  2B25aP 

533  =  455^^  +  AuP^  +  cf  [Hna^  +  2{Hi2  +  2Hm)a‘^p'^  +  H22p^] 

534  =  -^550!  —  Cl  +  (A2  +  2.F66)n:)S^j 

535  =  -<4440:  —  Cl  F22P^  +  (A2  +  2.F66)c:^/3j 

S44  =  .A55  +  .Diia^  +  Deep^,  S45  =  (5i2  +  DeQ)aP,  S55  =  .A44  +  Dqqo?  +  jD22/3^710) 
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The  mass  coefficients  are  the  same  as  those  defined  in  Eq.  (7.6b),  except  for  77134  =  77735  =  0. 
The  vector  {E}  is  given  by 


0 

0 

mLn 

l3NLn 

^  Qmn 

^  —  < 

0 

>  -  < 

0  ► 

0 

^•Vlmn 

«^lmn 

0  . 

'•  P^2mn  ‘ 

^  P^2mn  > 

(7.11) 


8.  Finite  element  formulation  of  CLPT 


8.1.  Weak  forms 


In  this  section,  finite  element  models  of  Eqs.(3.10)-(3.12)  governing  the  motion  of 
laminated  plates  with  piezoelectric  layers  are  developed  for  the  linear  case  first.  Then  the 
nonlinear  finite  element  model  will  be  discussed.  The  finite  element  model  is  based  on  the 
weak  forms  of  Eqs.(3.10)-(3.12).  The  weak  forms  are  (see  Reddy^^,  Chapter  10) 


^  f  \d6uQ^^  ,  ,  re  ^^“0  re  ^  fdwo\ 

°  “  /ne  [  dx  +  dy  dt^  (  ax  j 

-  {Nxxnx  +  NxyUy)  Sulq  ds 

0  =  1 


dxdy 


(8.1a) 


[  dSvo  a^vo 


[  dx 


^  r  e  ^^^0  r  J  fdwo\ 

dy  ~  ) 


dxdy 


(8.1!,) 


=  /  f- 


'0  d'^Swo  - . 

■^xx  2  ^  ^  j\dxv  r\  n  ■^yy  uW{)Q 


d‘^8w\ 

dx^  dxdy  dy"^ 

d^wo  ,  -  f  dSwo  d^WQ  dSwQ  d^WQ  \ 

+ + h 

/  d6wo  d^UQ  dSwQ  a^vo  \ 


I - v  , - -  -V  I 

dx  dt^  dy  dt"^  j 


dxdy  —  ^  VnSwQ  ds 


+ 


/  r  dSwo  . , ,  , ,  .  d6wo  . . ,  ,  ^  ' 

^  (Afxx^x  “k  j^xy'^y)  “k  ^  (jklxy^x  “k  ^^yy'^y) 


ds  (8.1c) 


where  {nx,ny)  denote  the  direction  cosines  of  the  unit  normal  on  the  element  boundary  F® 
of  the  element,  (^uoj^’yoj^iyo)  are  the  virtual  displacements,  which  take  the  role  of  weight 
functions,  and  14  is  the  effective  shear  force  defined  in  Eq.  (2.16).  The  stress  and  moment 
resultants  Nxx,Mxx:  etc.  are  known  in  terms  of  the  displacements  {v{),vo,wo)  and  thermal 
and  piezoelectric  forces  and  moments  through  Eqs.  (3.30)-(3.33). 


8.2.  Spatial  approximations 

First,  we  note  that  the  stress  and  moment  resultants  contain  first-order  derivatives 
of  (uq,vo)  and  second-order  derivatives  of  wq  with  respect  to  the  coordinates  x  and  y. 
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Second,  the  primary  variables  'Uq?  dw^/dx^  and  dw^/dy  must  be  carried  as  the  nodal 

variables  in  order  to  enforce  their  interelement  continuity.  Thus,  the  displacements  {uq^vq) 
must  be  approximated  using  the  Lagrange  interpolation  functions,  whereas  wq  should  be 
approximated  using  Hermite  interpolation  functions  over  an  element  fl^.  Let 

m 

uo{x,y,t)  « 

j=l 
m 

VQ{x,y,t)  w 

j=i 
n 

WQ{x,y,t)  «  A|(iV|(x,y) 

/C=l 

where  {u^Vj)  denote  the  values  of  (txoj'^o)  at  the  jth  node  of  the  Lagrange  elements,  A| 
denote  the  values  of  wq  and  its  derivatives  with  respect  to  x  and  y  at  the  fcth  node,  and 
Lagrange  and  Hermite  interpolation  functions,  respectively. 

There  exists  vast  literature  on  triangular  and  rectangular  plate  bending  finite  elements 
of  isotropic  or  orthotropic  plates  based  on  the  classical  plate  theory^^.  Here  we  discuss 
rectangular  plate  bending  elements.  There  are  two  kinds  of  plate  bending  elements. 
A  conforming  element  is  one  in  which  the  interelement  continuity  of  9x  =  dwo/dx^  and 
9y  =  dwo/dy  (or  dwo/dn)  are  satisfied,  and  a  nonconforming  element  is  one  in  which  the 
continuity  of  the  normal  slope,  dw^/dn^  is  not  satisfied  (see  Figure  8.1). 


(8.2a) 

(8.26) 

(8.2c) 


Figure  8.1:  (a)  Nonconforming  and  (b)  conforming  finite  elements  for  plate  bending 
based  the  classical  plate  theory. 
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A  nonconforming  rectangular  element  has  wq,  9x^  and  9y  as  the  nodal  variables  (see 
Figure  8.1a).  The  normal  slope  varies  cubically  along  an  edge  whereas  there  are  only  two 
values  of  dwo/dn  available  on  the  edge.  Therefore,  the  cubic  polynomial  for  the  normal 
derivative  of  wq  is  not  the  same  on  the  edge  common  to  two  elements.  The  interpolation 
functions  for  this  element  can  be  expressed  compactly  as 

V’i  (*  =  1,4,7,10);  (pf  =  gi2  (i  =  2,5,8,11)-,  v?|  =  5*3  (i  =  3, 6, 9, 12)  (8.3a) 

where 

9ii  =  ^(1  +  Co)(l  +  %)(2  +  io  +  Vo-f-  rf) 

gi2  =  g4i(^0  —  1)(1  +  %)(1  +  ^0)^,  9i3  =  “  1)(1  +  ^o)(l  +  Vo)^ 

^^{x-  Xc)/2,  •n={y-  yc)/b,  Co  =  CCu  Vo  =  WH  (8.36) 

where  2a  and  26  are  the  sides  of  the  rectangle,  and  {xc,yc)  are  the  global  coordinates  of  the 
center  of  the  rectangle. 

A  conforming  rectangular  element  with  wq,  dwo/dx,  dwo/dy,  and  d'^wo/dxdy  as  the 
nodal  variables  is  shown  in  Figure  8.1b.  The  interpolation  functions  for  this  element  are 

Pt  =  9n  {i  =  l,  5, 9, 13);  pf  =  9i2  {i  =  2, 6, 10, 14) 

Pi  =  9i3  {i  =  3, 7, 11, 15);  p^  =  ^i4  {i  =  4, 8, 12, 16)  (8.4o) 

where 

9il  “  Ci)^(Co  -  2)(tJ  +  T}i)^(T}0  -  2),  9i2  =  +  6)^(1  ~  Co)(^  +  Vi)‘^{V0  -  2) 

9i3  =  ^ViiC  +  Ci)^(Co  -  2)(77  +  r?i)2(l  -  rio),  +  Ci)^{l  -  Cq){v  +  vd'^il  “  ^o) 

(8.46) 

The  conforming  element  has  a  total  of  six  degrees  of  freedom  per  node,  whereas  the 
nonconforming  element  has  a  total  of  five  degrees  of  fireedom  per  node.  For  the  conforming 
rectangular  element  (m  =  4  and  n  =  12)  the  total  number  of  nodal  degrees  of  freedom  per 
element  is  24,  and  the  nonconforming  element  the  total  number  of  degrees  of  freedom  per 
element  is  12. 

8.3.  Semidiscrete  finite  element  model 

Substituting  approximations  (8.2)  for  the  displacements  and  the  ith  interpolation 
function  for  the  virtual  displacement  {6xiq  ~  ‘ipi,6vo  ~  rpi,6wo  ~  pi)  into  the  weak  forms,  we 
obtain  the  ith  equation  associated  with  each  weak  form 


m 

o  =  E 

+  K}fv]  +  Mnufj  + 

rpl  rpTl 

^  i  i 

-fP 

j=l 

k=l 

m 

n 

o  =  E 

+  M^i)  +  ^  (Kf  A|  +  Mf  A|) 

—  Ff  — 

-fP 

3=1 

k=i 

m 

n 

o-E 

(Kgu|  +  Kl]v]  +  Mgu^j  +  M§v^}  +  A|  +  M|| A|) 

3=1 

e=i 

-fI-fP-fP 

(8.5) 
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where  i  =  1, 2,  •  •  • ,  m;  A;  =  1, 2,  •  •  • ,  n.  The  coefficients  of  the  stiffness  matrix 

mass  matrix  (symmetric),  and  force  vectors  F“,  and  Ff°‘  are  defined  as 

follows: 

K}}  =  AnSf?  +  ^16  {sfy  +  5r;)  +  Aeesyy 

^ij  —  +  ^165^^  +  A2eSfj  +  Aqq^j 

K‘if  =  JkeSf^  +  ^26  {Sff  +  5^")  +  A22S!]' 

kE  =  -Sniir  -  -  2B2eRr  "  -  B2,RT  "  2^66/?^ 

Ki  =  -BnRT  -  B26l^r  -  2BeeR-k^  -  -  B22RT  "  2B26Rr 

Klj  =  +  Di2  +  2Di6 

+ 2i>26  {Ttr^ + Ttr^) + 4^667^r^ + D22itr^ 

M]f  =  lo^trj  dxdy,  ^ 

=  jf  ^  /oV-l^l  dxdy,  Mi  =  -  jT  ^  dxdy 

—  (p  ij^xx'^x  “1“  ^xy'^y)  ds^  Fj^  —  I  (^^xy'^x  “1“  ^yy'^y^  ds 

Jr^  JTe 

n  =  j^  v%  dxdy  +  £  (v„v% + r.M  +  r,M)  * 


dx  dv 


g^2  -t-  -t-  Qy2 


= L  ^  = L  ( fa  "'S  +  % 

=  L  +  ^K)  dxdy 

= L  "  f  = £  (f  <. + f 

=  fa. 


Br= 


dxdy  dy^  yy 


(8.6a) 


r  dr.d-^ifi 


(8.66) 


hed^dvdCdfi — " 

and  ^,»7,C)  and  /x  can  be  equal  to  x  or  y.  In  matrix  notation,  Eq.  (8.5)  can  be  expressed  as 
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8.5.  Nonlinear  finite  element  model 

The  nonlinear  finite  element  model  has  the  form 

[K^]{A^}  +  [M^]{A}  =  {F^}  (8.8a) 

[K^]  =  [Kl  +  +  K%r  +  K%p]  (8.86) 

where  [Kf]  is  the  linear  stiffness  matrix,  and  [Kffj]  is  the  geometric  nonlinear  part,  [Kffj] 
the  thermal  nonlinear  part,  and  [Kffp]  the  piezoelectric  nonlinear  part  of  the  stiffness  matrix 
due  to  the  nonlinear  terms  in  the  equations  of  motion.  The  elements  of  the  linear  stiffness 
matrix  [Kf]  and  mass  matrix  [M®],  which  is  linear,  were  given  in  Eqs.(8.6a,b).  The  nonzero 
nonlinear  stiffness  coefficients  \K°^\  (a  =  1,2, 3)  are 


i^N  .  • 


-\L 

'  dy 


N6j  +  ^N2j^  dxdy 
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\  dx  dy 


13 

ji 
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(dwp's^ 
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fdw^ 


Ula{^)  +A2a(^)  +2^„ 
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\dx  J 
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(8.9o) 


(8.96) 


y^NT^ 

1 

{k^p 

)..=  / 

\  / 

JOfi  1 

1 1  -  -S  (t  t t  ^)] 


for  a  =  1,2,6.  The  same  expressions  hold  for  the  first-order  plate  theory  with  (pi  replaced 
by  il)i,  because  wq  is  also  approximated  by  the  Lagrange  functions  ipi  in  FSDT.  We  note  that 
the  nonlinearity  of  the  coefficients  and  is  solely  due  to  the 

transverse  deflection,  too,  as  can  be  seen  from  Eqs.  (8.9a,  6).  In  addition,  the  stiflEness  matrix 
[K^]  is  not  symmetric  for  the  nonlinear  case  because 

Kl^  =  2Kf,  a  =  1,2  (8.12) 


The  nonlinear  algebraic  equations  of  the  bending  problem  or  fully  discretized  transient 
problem  must  be  solved  by  an  iterative  method.  In  iterative  methods,  the  nonlinear  equations 
are  linearized  by  evaluating  the  nonlinear  terms  with  the  known  solution  from  preceding 
iteration(s).  Two  commonly  used  iterative  methods  are:  (1)  the  Picard  method,  and  (2)  the 
Newton-Raphson  method.  While  the  Picard  method  is  simple  both  in  concept  and  computer 
implementation,  the  method  uses  the  stiffness  matrix  [k%  which  is  unsymmetric.  The 
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Newton-Raphson  method  is  based  on  the  Taylor  series  expansion,  and  it  uses  the  tangent 
stiffness  matrix,  which  is  symmetric  for  all  structural  problems.  Also,  the  method  has  faster 
convergence  for  most  applications  than  the  Picard  method.  Here  a  brief  discussion  of  these 
two  iterative  methods  is  presented.  For  the  sake  of  discussion,  we  use  the  fully  discretized 
equation 

=  (8.13) 

is  used. 

In  the  Picard  method,  also  known  as  the  direct  iteration  method,  the  solution  vector 
from  the  previous  iteration  is  used  to  evaluate  the  stiffness  matrix,  and  the  solution  at  the 
subsequent  iteration  is  determined  by  solving  the  assembled  equations  after  the  imposition 
of  boundary  and  initial  conditions.  At  the  element  level,  the  Picard  iteration  scheme  may  be 
expressed  as 

=  {F^}  (8.14) 

where  {A®}^  denotes  the  solution  vector  at  the  rth  iteration.  Thus,  in  the  direct  iteration 
method,  the  coefficients  Kfj  are  obtained  by  evaluating  them  with  wo{x,  y,  t)  from  the  rth 
iteration.  At  the  beginning  of  the  iteration  {i.e.,  r  =  0),  we  assume  that  {A®}°  =  {0}  so  that 
the  solution  at  the  first  iteration  is  the  linear  solution,  because  the  nonlinear  stiffness  matrix 
reduces  to  the  linear  one.  The  iteration  process  is  continued  until  the  difference  between 
{A®}’'  and  {A®}’’"^^  reduces  to  a  preselected  error  tolerance.  The  global  error  criterion  is  of 
the  form 


<  e  (say  10  (8.15) 


where  N  is  the  total  number  of  nodal  generalized  displacements  in  the  finite  element  mesh, 
and  e  is  the  error  tolerance.  All  quantities  in  Eq.  (8.15)  are  understood  to  be  at  the  global 
level. 

The  Newton-Raphson  iterative  method  is  based  on  Taylor’s  series  expansion  of  the 
nonlinear  algebraic  equation  (8.13)  about  the  known  solution.  Suppose  that  the  solution 
at  the  rth  iteration,  {A®}’',  is  known.  Let 

{R®}  =  [.^®]{A®}-{.F®}  =  0  (8.16) 

where  {R®}  is  called  the  residual,  which  is  a  nonlinear  function  of  the  unknown  solution  {A®}. 
Expanding  {i?®}  in  Taylor’s  series  about  {A®}'',  we  obtain 


{0}  =  {B*}  =  {IfY  +  ({A‘r+'  -  {A'D 


1  r  a2{R®} 


2!  a{A«}a{A^} 


({A®}»-+1  -  {A®}’’)^  +  ... 


(8.17o) 


or 


0  =  {R®}’-  +  {5A®}  +  0({6A®}2) 


(8.176) 
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where  O(-)  denotes  the  higher-order  terms  in  and  is  the  tangent  stiffness 

matrix  (or  geometric  stiffness  matrix) 


/  ^  \  tan 

{ \K%]  = 

d{R^y 

jr  /  — 

evaluated  at  {A®}  = 


and  the  increment  of  the  solution  vector  is  defined  by 

{^A^}  =  {A"}^+i  -  {A^Y 


(8.18) 


(8.19) 


Equation  (8.17b)  is  approximated  by  neglecting  terms  of  order  2  and  higher  in  the  solution 
increment  {^A®}.  We  obtain 


([i:^({A^}n])*“”  {^A^}  =  (8.20) 

The  assembled  equations  are  then  solved  after  imposing  the  boundary  and  initial  conditions 
of  the  problem.  We  have 

{«A}  =  -  {flr 

=  (|A-({A}'-)1“”)'‘  ({#}  -  [A:({A}--)l{Ar)  (8.21) 

and  the  total  solution  at  the  (r  +  l)th  iteration  is  given  by 

{^}r+l  ^  (g  22) 

The  iteration  process  is  continued  by  solving  Eq.  (8.20)  until  the  convergence  criteria  in  Eq. 

(8.15)  is  satisfied  or  the  residual  {R}j  measured  in  the  same  way  as  the  solution  error  in  Eq. 

(8.15) ,  satisfies  the  error  criterion.  Note  that  at  the  beginning  of  each  iteration  the  tangent 
stiffness  matrix  and  residual  vector  must  be  updated  using  the  latest  available  solution  {A}. 
If  the  tangent  stiffness  matrix  is  kept  constant  for  a  preselected  number  of  iterations  but  the 
residual  vector  is  updated  during  the  iteration,  there  can  be  a  computational  saving.  Such 
an  approach  is  called  the  modified  Newton-Raphson  method. 

9.  Finite  element  model  of  TSDT 
9.1.  Introduction 

The  primary  variables  of  the  third-order  theory  are  {un^UsjWo^WQ^n  =  d'^o/dn,(j)n^(j)s)j 
where  (un^Us)  denote  in-plane  normal  and  tangential  displacements,  and  ((^nj^s)  are  the 
rotations  of  a  transverse  line  about  the  in-plane  normal  and  tangent.  A  displacement 
finite  element  model  based  on  Eqs.  (4.12)-(4.16)  requires  the  Lagrange  interpolation  of 
{uQ^voy  <j>x^  <l>y)  and  Hermite  interpolation  of  tUQ.  A  conforming  element  will  have  eight  degrees 
of  freedom  whereas  a  nonconforming  element  will  have 

Seven  degrees  of  freedom  per  node.  In  view  of  the  detailed 
discussion  of  finite  element  models  of  the  classical  plate  theory  presented  in  Section  8,  only 
the  salient  features  of  the  model  are  discussed  here^^. 
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9.2.  Finite  Element  Model 

The  generalized  displacements  are  approximated  over  an  element  by  the  expressions 

m 

uo{x,y,t)  =  Y^uf{t)'ipf{x,y) 

i—l 

m 

vo{x,y,t)  =  ^  (<)■!/'•  (a;,  y) 

i=l 

m 

wo{x, y,t)  =  ^Af(t)(/?f(a;,y) 
i=l 
m 

<l>x{x,y,t)  =  Y^Xt{t)'ipl{x,y) 

(9.1) 


1=1 

m 


4>y{x,y,t)  =  Y^Yt{t)ipi{x,y) 
i=l 

where  tpf  denote  the  Lagrange  interpolation  functions  and  are  the  Hermite  interpolation 
functions.  Here  we  chose  the  same  approximation  for  the  in-plane  displacements  {uq,vo)  and 
rotations  {(f>x,(l>y),  although  one  could  use  different  approximations  for  these  two  pairs.  In 
the  case  of  the  conforming  element,  the  four  nodal  values  associated  with  wq  are 


Ai  =  too,  A2  = 


dwo 

dx 


A3  = 


dwQ 

dy 


A4  = 


5^  too 
dxdy 


(9.2) 


For  the  nonconforming  element,  the  cross  derivative  is  omitted.  The  conforming  rectangular 
element  with  linear  interpolation  of  the  in-plane  displacements  and  rotations  has  eight  degrees 
of  freedom  per  node.  The  corresponding  nonconforming  element  has  seven  degrees  of  freedom 
per  node. 

Substitution  of  Eq.  (9.1)  into  the  weak  form  of  Eqs.  (3.12)-(3.16)  yields  the  finite  element 
model 
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or,  in  compact  form,  we  can  write 
5 

=  i  =  l,2,...,na  (9.36) 

/3=1  j=l 

where  a  =  1, 2, 3, 4, 5;  ni  =  712  =  n4  =  ns  =  4  and  —  16  for  the  conforming  element.  The 
nodal  values  Af ,  the  linear  stiffness  coefficients  Kff,  and  mass  coefficients  Mff  are  defined 
by 

A]  =  vj,  A?  =  Ai,  A|=X',  Af  =  i^?  (9.4) 
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where  Vn  is  defined  in  Eq.  (3.21),  and  P„„  is  related  to  Pxx,Pxy,  and  Pyy  in  the  same  way 
Mnn  related  to  Mxx,Mxy,  and  Myy  [see  Eq.  (2.186)].  The  following  notation  is  used  in 
defining  the  above  coeificients: 

cO  f  oOl  _  f  j„j.. 
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(9.9) 


(9.10) 


The  nonlinear  stiffness  matrices  remain  the  same  as  given  in  Eqs.  (9.9)-(9.11). 

Note  that  the  stiffness  matrix  evaluation  requires  the  computation  of  the  second 
derivatives  of  the  interpolation  functions  used  for  the  transverse  deflection.  The  isoparametric 
elements  are  used  for  {uo,vo,<px,^y)  and  subparametric  formulation  is  used  for  wq.  The 
transformation  equations  required  to  numerically  evaluate  the  stiffness  coefficients  can  be 
found  in  the  textbooks  by  Reddy®2>53  xhe  discussion  presented  there  on  shear  locking 
with  the  shear  deformable  displacement  finite  element  model  also  applies  to  the  present 
element  in  evaluating  the  shear  stiflhiess  coefficients,  which  should  be  evaluated  using  reduced 
integration. 


In  the  case  of  the  first-order  shear  deformation  plate  theory,  the  resulting  finite  element 
models  require  only  C®— continuity  of  all  generalized  displacements  {uo,vo,WQ,(l>x,(l)y),  and 
they  are  the  most  economical  while  accounting  for  the  transverse  shear  deformation.  The 
theory  is  discussed  in  the  next  section.  However,  the  finite  element  models  will  not  be 
discussed^^ 
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10.  Analysis  of  magnetostrictive  plates  using  FSDT 
10.1.  Introduction 

The  grains  of  certain  materials  consist  of  numerous  small,  randomly  oriented  magnetic 
domains  that  can  rotate  and  align  under  thew  influence  of  an  external  magnetic  field.  The 
magnetic  orientation  brings  about  internal  strains  in  the  material.  This  is  known  as  the 
magnetostriction.  This  magneto-elastic  response  is  non-linear.  A  commercially  available 
magnetostrictive  material  Terfenol-D  is  an  alloy  of  terbium,  iron,  and  dysprosium. 

Laminated  composite  plates  containing  magnetostrictive  layers  are  modelled  as 
distributed  parameter  systems  and  the  magnetostrictive  layers  are  used  to  control  the 
vibration  suppression^®"®®.  Velocity  feedback  with  constant  gain  distributed  controller  is 
chosen  to  achieve  vibration  suppression. 


10.2  Governing  equations 

Consider  a  symmetrically  laminated  composite  plate  of  n  layers  with  the  mth  and 
(n  —  m  -|-  l)th  layers  being  made  of  magnetostrictive  material,  and  the  remaining  n  —  2  layers 
made  of  a  fiber-reinforced  material  and  having  varying  fiber  orientation  6  and  symmetrically 
disposed  about  the  center  plane  of  the  plate  (see  Figure  10.1).  The  displacement  field  of  the 
first-order  shear  deformation  theory  (FSDT)  is  of  the  form 


u{x,y,z,t)  =  uo{x,y,t)  +  z^x{x,y,t) 
v{x,y,z,t)  =  vo{x,y,t)  +  z<f)y{x,y,t) 

w{x,y,z,t)  =  wo{x, y,t)  (10.1) 


where  {uo,vo,wo,(f>x,(j)y)  are  unknown  functions  to  be  determined,  {v.o,vo,wo)  denote  the 
displacements  of  a  point  on  the  plane  z  =  0,  and  which  indicate  that  (f>x  and  <j)y  are  the 
rotations  of  a  transverse  normal  about  the  y—  and  x— axes,  respectively. 
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A  «th  Layer  (rt-/n+l)th  (magnetostrictive)Layer 


'l^4s 

1X1 

"  <jr'7!l'77:7 

1st  Layer  mth  (magnetostrictive)Layer 


Figure  10.1:  Symmetrically  laminated  plate  with  magnetostrictive  layers. 
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The  linear  strains  associated  with  the  displacement  field  (10.1)  are 
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(10.2) 


The  constitutive  relations  of  the  /cth  lamina  are 
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(10.3a) 

(10.3b) 


where  are  the  plane  stress-reduced  stiffnesses  and  are  the  magneto-mechanical 
coupling  moduli  of  the  fcth  lamina  [see  Eqs.  (3.6)-(3.8)  and  (4.8)]. 

The  governing  equations  of  the  first-order  theory  can  be  derived  using  the  dynamic 
version  of  the  principle  of  virtual  displacements.  The  equations  of  motion  are 


dNxx  ,  dNxy  _  J  ,  r 


dNxy  dNyy _  d'^VQ  d‘^<j)y 

dx  dy  “  °  at2 

dQx  ,  dQy  d^wo 


0]\^xx  , 

dx  dy 

d]\/lxy  ^^^yy 
dx  dy 


Qx  — 12 


d^(t>, 


X  ,  J  df^UQ 

+  h' 


dt^ 

-Qy  =  -^2-^  +  h 


df^ 

d^vo 

df^ 


The  natural  boundary  conditions  are  of  the  form 


^nn  ^nn  —  0  5  ^ns  ^ns  —  0  ,  Qn  Qn  0 
^nn  ^nn  “  0  >  ^ns  ^ns  ~  0 


where 


Qn  =  Qx'^x  “h  Qy'^y 


Thus  the  primary  and  secondary  variables  of  the  theory  are 


primary  variables:  wq^  <j)n^  (f)s 

secondary  variables:  Nnny  Nns^  Qny  Mnn^  Mns 


(10.4a) 

(10.4b) 

(10.4c) 

(10.4d) 

(10.4e) 


(10.5) 

(10.6) 

(10.7) 


39 


The  initial  conditions  of  the  theory  involve  specifying  the  values  of  the  displacements 
and  velocities  at  time  t  =  0: 


Un  =U^,  Us  =  Wo  =  W^,  (pn  =  0°,  <l>s  = 

iln  ^Un,  Us  =  Ws,  WO  =  Wq,  K  =  <^n.  (10-8) 


for  all  points  in  Qq- 


10.3  Velocity  feedback  control 

Considering  velocity  proportional  closed-loop  feedback  control,  the  magnetic  field 
intensity  Hz  is  expressed  in  terms  of  coil  current  I  {x,  y,  t)  as 

H{x,y,t)  =  kclix,t)  (10.9) 

and  I{x,y,t)  is  related  the  velocity  wq  by 


Iix,y,t)  =  c{t)^  (10.10) 

where  kc  is  the  coil  constant,  which  can  be  expressed  in  terms  of  the  coil  width  coil  radius 
To,  and  number  of  turns  ric  in  the  coil  by 


kc  — 


Uc 

+  4r^ 


(10.11) 


and  c{t)  is  the  control  gain. 

10.4  Laminate  constitutive  equations 

The  force  and  moment  resultants  are  related  to  the  strains  by 


where  K  is  the  shear  correction  factor  and 

(^^fc+1  -  Zk) 

Jfc=l 

Dij  =  \j2Qf  (4+1-4) 

(A44,  ^45,^455)  =  y  (Q^^ ,  055^)('2^*:+1  ~  ^k) 

k=l 


(10.12a) 

(10.12b) 


(10.13) 
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and  the  magnetostrictive  stress  resultants  {N^}  and  {M^}  are  defined  by 

{S}=  s 

yy  ^  k=m,n-m+l  *•  J 


(10.14) 


zHz  dz 


=  ckc  E  f" 

fc=m,n— m+l 


dwo 

dt 


z  dz 


f  Esi  1  dwQ 

1^32  J  dt 


(10.15) 


10.5  Anedytical  solution 

Here  we  seek  the  Navier  solutions  for  the  case  of  simply  supported,  symmetric  laminates. 
Towards  this  end,  we  write  the  governing  equations  in  terms  of  the  generalized  displacements 

{wq,  <f)x,  <i>y)'- 


dx 


_a 

dy 


d^wo 

df^ 


( 


OX 


+  D^f-  (^  +  ^) 

OX  \  oy  ox  J 

d^(t>x 


\  c  ^^0 


dt^ 


d<px 

dy 


dy 


(10.16) 


(10.17) 


(10.18) 


The  simply  supported  boundary  conditions  for  the  first-order  shear  deformation  plate 
theory  (FSDT)  are 

<f>x{x,0,t)  =  0,  <l)x{x,b,t)  =  0,  ^y{0,y,t)  =  0,  ,  (l>y{a,y,t)  =  0 

wo(x,0,t)  =  0,  wo{x, b,t)  =  0,  u;o(0,y,t)  =  0,  wo{a,y,t)  =  0 
Mxx{0,y,t)  =  0,  Mxx{a,y,t)  =  0,  Myy(x,0,t)  =  0,  Myy{x,b,t)  =  0  (10.19) 
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The  boundary  conditions  in  (10.19)  are  satisfied  by  the  following  expansions 

OO  00 

wo{x,y,t)  =  53  Wmn{t)  sin  ax  sin  Py  (10.20a) 

n=l  m=l 

OO  OO 

Px{x,y,t)  =  53  53  Xmn{t)  cos  ax  sin  Py  (10.20b) 

n—1  m—1 

OO  OO 

<f)y  {x,  y,  0  =  53  53  ^rnnit)  sin  ax  cos  Py  (10.20c) 

n=l  m=l 


The  mechanical  load  and  magnetostrictive  moments  are  also  expanded  in  double  Fourier 
sine  series 

OO  OO 

q{x,  y,  0  =  53  53  Qmn{t)  sin  ax  sin  py  (10.21a) 

72=1  772=1 
OO  OO 

M^{x, y, *)  =  53  53  sin axsinPy  (10.21b) 

72=1  772=1 
OO  OO 

{x,  y>  ^)  =  5Z  53  sin  ax  sin  py  (10.21c) 

n=l  772=1 


where,  for  example, 


/I  4:  r  rri'KX  .  nTry 

Qmn{t)  =  “T  /  /  sin -  sin-—  dxdy 

do  Jo  Jo  a  0 

=  /  -^^(^.y^^)  sin^^  sin^^  dxdy 


(10.22a) 

(10.22b) 


Substitution  of  Eqs.  (10.20)-(10.22)  into  Eqs.  (10.16)-(10.18)  yields  the  equations 

533  S34  %  rWmra'I  0  0  O'  f  IFmri  ^ 

534  S44  S45  \  ^mn  /  "I"  I  C12  0  I  )  ^mn  ? 

.S35  S45  S55.  V  ^mn  )  ,^13  0  0_  Yfnn  } 

'mss  0  0  1  f  ’ifmn  I  (  Qmn  ] 

+  0  iflu  0  <  Xmn  >  =  •{  0  ?  (10.23) 

.  0  0  rhssj  [  Ymn  )  I  0  J 

where  Sij,  Cij,  and  rhij  are  defined  by 

S33  =  K{AssoP'  +  A44P^) 

I34  =  KAssa,  535  =  KAuP,  S44  =  {Dncn^  +  DqsP^  +  KA^s) 

S45  =  (Di2  +  Des)aP,  S55  —  {Deea^  +  Z?22/?^  +  KA^/p 
C\2  =  Q!^12,  Ci3  =  P£\s 
msi  =  Iq,  mu  =  h,  =  h 


(10.24a) 

(10.24b) 


where  the  magnetostrictive  coefficients  5i2,  and  £13  are  defined  in  Eq.  (10.18). 
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For  vibration  control,  we  assume  g  =  0  and  solution  of  the  ordinary  differential  equations 
in  Eq.  (10.26)  in  the  form 

Wmnit)=Woe^\  Xmnit)=Xoe^\  Ymn{t)  (10.25) 

and  obtain,  for  non-trivial  solution,  the  result 

■§33  S34,  S35 

543  ^44  1S45  =  0  (10.26) 

553  <554  <555 

where 

Sij  =  Sij  +  XCij  +  Mij  (10.27) 

where  i,j  =  3,4,5.  This  equation  gives  three  sets  of  eigenvalues.  The  lowest  one  corresponds 
to  the  transverse  motion.  The  eigenvalue  can  be  written  as  A  =  —a+iuJa,  so  that  the  damped 
motion  is  given  by 

/  1  _Q,t  .  ,  .  ^ 

wo{x,y,t)  =  — e  sinwdisin - sin— ^  (10.28) 

a  u 

In  arriving  at  the  last  solution,  the  following  initial  conditions  are  used; 

'u;o(2;,2/,0)  =  0,  wo{x,y,0)  =  l,  ^&i(a:,y,0)  =  0,  ^x{x,y,0)  =  0 

<j)y{x,y,0)  =  0,  ^y{x,y,0)  =  0  (10.29) 

10.6.  Numerical  results  and  discussion 

Numerical  studies  were  carried  out  to  analyze  damped  natural  frequencies,  damping 
coefficients,  and  the  vibration  suppression  time,  using  the  three  theories.  Different  lay-ups 
were  used  to  show  the  influence  of  the  position  of  magnetostrictive  layer  on  the  vibration 
suppression  time.  A  time  ratio  relation  between  the  thickness  of  the  layers  and  the  distance 
to  the  neutral  axis  of  the  laminated  composite  beam  is  found.  All  values  of  the  material  and 
structural  constants  are  indicated  in  the  tables. 

The  numerical  values  of  various  coefficients  (namely,  the  inertial  and  magnetostrictive 
coefficients)  based  on  different  lay-ups  and  material  properties  [CFRP,  Graphite-Epoxy  (AS), 
Glass-Epoxy  and  Boron-Epoxy]  are  listed  in  Tables  10.1  and  10.2.  Table  10.2  also  shows  the 
damping  coefficients  and  natural  frequencies  for  different  materials  and  lay-ups.  The  damping 
and  frequency  parameters  for  transverse  modes  n  =  1  to  n  =  5  are  shown  in  Table  10.3,  and 
they  are  compared  with  the  results  obtained  by  Murthy  et  al.  [56]  using  the  Euler-Bernoulli 
beam  theory  (EBT).  Only  in  the  higher  modes  there  is  some  difference  between  the  numerical 
results  predicted  by  the  three  theories.  Table  10.4  shows  the  influence  of  the  position  of  the 
magnetostrictive  layer  in  the  2?-direction  and  the  influence  of  the  lamination  scheme  in  the 
damping  and  frequency  parameters.  The  value  of  a  increases  when  the  magnetostrictive  layer 
is  located  further  away  from  the  ar-axis,  indicating  faster  vibration  suppression.  The  lay-up 
[m/904]5  represents  the  softest  beam  and  the  lay-up  [m/04]s  the  stiffest  beam. 
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Table  10.1.  Coefficients  for  different  Lay-ups  and  materials. 


Material 

Lay-up 

D„  (10’) 

H„  (10-’) 

AssOO*) 

FssdO-’) 

3.739 

9.333 

6.620 

6.902 

[45/in/-45/0/90]s 

3.552 

[IKEZI 

6.620 

HBEI 

8.792 

CFRP 

[m/±45/0/90]s 

3.303 

4.069 

13.168 

[in/904]s 

1.432 

2.567 

5.063 

wmmt 

13.168 

[m/OJs 

7.015 

7.927 

11.189 

6.620 

7.506 

13.168 

Gr.-Ep  (AS) 

[±45/m/0/90]s 

5.629 

beesi 

Gl.-EP 

[+45/tn/0/90]s 

3.700 

■IQQ 

Br.-Ep 

[±45/ni/0/90]s 

5.730 

8.259 

14.865 

5.634 

7.569 

Data:  CFRP :  Eii=138.6Gpa  E22=8.27  Gpa  Gi3=G23=0.6  E22  Gi2=4.12Gpa  Vi2=0.26  p=1824kg.m‘^ 

Graphite-Epoxy  (AS)  : El  1=1 37.9  GPa  £22=8.96  Gpa  Gi2=Gi3=7.10  Gpa  G23=6.21  Gpa  Vi2=0.30  p=1450kg.m'^ 

Glass-Epoxy:  Eii=53.78  Gpa  E22=17.93  Gpa  Gi2=Gi3=8.96  Gpa  G23=3.45  Gpa  Vi2=0.25  p=1900kg.m'^ 

Boron-Epoxy:  En=:206.9Gpa  £22=20.69  Gpa  Gi2=  Gi3=6.9  Gpa  G23=4.14Gpa  Vi2=0.30  p=:1950kg.m‘^ 


Table  10.2.  Inertial  and  magnetostrictive  coefHcients  and  the  parameters  a  and 

Q)dii. _ 


Material 

Lay-up 

lo 

l2(10^) 

l4(10-’) 

l6(10'^) 

-e  (10"*) 

CFRP 

[+45/m/0/90]s 

33.092 

2.461 

2.907 

4.508 

22.128 

1.438 

3.30+104.85 

[45/m/-45/0/90]s 

33.092 

3.352 

4.600 

30.979 

3.872 

4.62+102.15 

[in/±45/0/90]s 

m*m]m 

8.521 

17.171 

8.165 

5.94+98.42 

[ni/904]s 

laaign 

■mi 

mm 

wm 

5.94±64.65 

[in/04]s 

33.092 

4.540 

mm 

■mi 

mm 

5.94+143.57 

Gr.-Ep 

[+45/ni/0/90]s 

Mugiril 

■mi 

mmm 

3.63±1 13.06 

Gl.-EP 

[±45/m/0/90]s 

mi 

■mi 

1.438 

3.24±85.54 

Br.-Ep 

[±45/in/0/90]s 

34.100 

4.782 

22.128 

1.438 

3.20+127.90 

Table  10.3.  Damping  and  frequency  parameters  due  to  the  transverse  modes. 


-OC  ±  COdn  (rad/s)  -  Lay-up  [±45/in/0/90]s 

Mode 

Murty  et  al 

EBT 

TBT 

RBT 

1 

3.29±104.88 

3.30+104.85 

3.301104.82 

3.30+104.82 

2 

13.19±419.50 

13.20+419.37 

13.17+418.90 

13.161418.80 

3 

29.70±943.88 

29.681943.40 

29.531941.05 

29.481940.52 

4 

52.86±1678.83 

52.73+1676.72 

52.2711669.32 

52.1011667.68 

5 

82.59+2621.87 

82.3412619.02 

81.2212601.04 

80.8012597.09 

Data:  CFRP  :  Eii=l 38.6  Gpa  E22==8.27  GPa  Gi2=4.12GPa  Gi3=G23=0.6  E22 

Magnetostrictive  layer  :  Eni=26.5  GPa  pm=9250  kg.m'^  dk=1.67xlO‘WA 

Vi2=0.26  p=1824kg.m*^ 

c(t).Rc=10'*  Vro=0  a=l  m 

Table  10.4.  Damping  and  frequency  parameters  for  different  lay-ups. 


-a  ±  COdn  (rad/s)  -  mode  1 

Lay-up 

Murty  et  al 

EBT 

TBT 

RBT 

[45/iii/-45/0/90]s 

4.60+102.17 

4.621102.15 

4.621102.12 

4.621102.11 

[in/145/0/90].s 

5.90198.44 

5.94198.42 

5.94198.39 

5.93198.38 

[nj/904ls 

5.90164.65 

5.94164.65 

5.94164.64 

5.94164.64 

[ni/04]s 

5.901143.58 

5.941143.57 

5.931143.49 

5.931143.44 
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A  comparision  of  the  fundamental  transverse  and  axial  modes,  obtained  using  the 
Timoshenko  beam  theory  (TBT)  and  Reddy  third-order  beam  theory  (RBT)  show  (results 
are  not  included  here  due  to  space  limitations)  that  the  Reddy  third-order  beam  theory 
give  lower  values  for  natural  frequency.  The  comparision  of  the  uncontrolled  and  controlled 
motion  at  the  mid  point  of  the  beam  is  shown  in  Figures  10.2-10.5  for  the  first  mode.  These 
figures  show  that  the  vibration  suppression  time  decreases  when  the  distance  to  the  neutral 
axis  is  increased,  and  it  remains  nearly  the  same  in  the  laminates  with  different  stiffness. 


Figure  10.2:  Comparison  of  original  and  controlled 
motion  at  the  midpoint  of  the  beam 
([45/-45/m/0/90]s). 


Figure  10.4;  Comparison  of  original  and  controlled 
motion  at  the  midpoint  of  the  beam 
([m/45/-45/0/90]s). 


Figure  10.3:  Comparison  of  original  and  controlled 
motion  at  the  midpoint  of  the  beam 
([45/m/-45/0/90]s). 


Time  (s) 

Figure  10.5:  Comparison  of  original  and  controlled 
motion  at  the  midpoint  of  the  beam 
([m/904]^). 
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